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ABSTRACT 


Elements ■with even atomic number (Z) in the interval 50 ^ Z ^ 50 
have been resolved in the cosmic radiation using the Heavy Nuclei 
Experiment on the HEAO-3 satellite. Their relative abundances have 
been compared with the results expected irom pure r-process materiel, 
pure s-process material, and solar system material, both with and 
without a modification due to possible first ionization potential effects. 
Such effects may be the result of the preferential acceleration, and 
hence enhancement in the cosmic rays, of those elements having low- 
first ionization potentials. We find that our measurements are incon- 
sistent with pure r-proeess material at the greater than 98% confidence 
level whether or not the first ionization potential adjustments are made. 

In addition, we have compared our results with mixtures having 
varying ratios of pure r-process material to pure s-process material. We 
find that, if no first ionization potential effects are included, 


(r/ s)cK s_ _ Q 2 o 
(r/s)ss 


+0.18 
— Q.]4 


where CRS refers to the cosmic ray source and S3 refers to the solar 
system, consistent with having an almost pure s-process source. If the 
first ionization potential adjustments are applied 


(r/ s) CR s 
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consistent with a solar system mixture. 
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I. INTRODUCTION 


Galactic cosmic rays are highly energetic particles in space, one of 
the few samples of extra-solar material available for direct study. As 
such, knowledge of their elemental and isotopic composition should ulti- 
mately help us to understand the processes of nucleosynthesis responsi- 
ble for their formation and for the formation of the chemical elements 
in general. For many lower charge nuclei the elemental composition of 
the cosmic radiation has already been reasonably well determined. For 
a much lesser number the isotope ratios have also been determined and, 
in some cases, show significant differences from the composition of the 
solar system (see e.g. Mewaldt 1903 and references therein). However, 
for approximately 2/3 of the periodic table, those elements heavier than 
iron (i.e. atomic number, Z. greater than 26), relatively little detail is 
known about the elemental, and nothing at all about the isotopic, cosmic 
ray abundances. Although comprising only a small fraction of the 
cosmic ray flux, and of the solar system these elements are important 
because the processes believed responsible for their formation consti- 
tute a distinct class of events which can best be studied in this charge 
range. This thesis will discuss measurements of a limited, but important, 
region from 5 oSn to S8 Ce, where differences between the possible 
nucleosynthesis mechanisms are particularly evident. 

Figure 1.1 is a graphic presentation of the abundances of the chemi- 
cal elements in the solar system as compiled by Anders and Ebihara 
(1962). The vertical axis is logarithmic in order to display the full range 
of variation present. Although these are the solar abundances (derived 
in large part from meteorite studies), to a rough first order the cosmic 
ray elemental abundances follow a similar curve. Note the general 
decline in abundance as one proceeds up scale, some 4-5 orders of mag- 
nitude r rom hydrogen to iron (Z=26), with an even sharper decrease 


Figure 1.1 


The abundances of the chemical elements in the solar system 
compiled by Anders and Ebihara (1882). 




immediately above iron. The ultra-heavy component, Z 30, is approxi- 
mately 4 orders of magnitude less abundant than iron. Note also the 
moderate abundance increase in the charge 50-60 region which falls off 
again above charge 50. 

Theories of nucleosynthesis, which attempt to explain the features 
of the solar system ana "cosmic" abundances in terms of nuclear 
processes occurring in stars, have been quite successful in reproducing 
the general features of this curve. One of the pioneering papers in this 
field is that of Burbidge, Burbidge, Fowler and Hoyle (1957). The forma- 
tion of the elements below iron is mainly the result of charged particle 
interactions which also serve as a stellar energy source. However, the 
synthesis of the elements above iron cannot be explained using these 
processes. By the time one reaches iron, the coulomb barrier between 
the interacting charged nuclei is so large that the temperature required 
to provide a particle with enough energy to penetrate the barrier will 
also disrupt the product nucleus. In addition, iron lies at the maximum 
in the binding energy per nucleon curve making formation of higher 
charged particles energetically unfavorable. It was Burbidge ei al 
(1957) who first clearly explained the formation of the elements above 
iron by means of neutron capture processes, thereby avoiding the 
coulomb barrier problems. 

These processes start with a "seed" nucleus, the result of prior 
nucleosynthesis, and a source of neutrons. The seed nuclei go through a 
sequence of neutron captures and beta decays continuing until the 
onset of fission or exhaustion of the neutron source. The seed is usually 
assumed to be a member of the iron group, the nuclides in the vicinity of 
the maximum in the binding energy per nucleon curve, most commonly 
M Fe. The source of neutrons is still a matter of some debate with the 
reactions I3 C(a,n) 16 0 and Z2 Ne(a,n) 25 Mg two likely candidates proposed 
by Cameron (1955, I960). The two main processes are actually the two 
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extreme eases o! a more general neutron capture process (Blake and 
Schramm 1976). These are the r (rapid) and the s (slow) process. 

In the s-process, the flux of neutrons is low enough that the time 
between subsequent neu'-ron captures is large compared with the beta 
decay time of any unstable nucleus encountered. As a result, the s- 
process follows a path along the valley of beta stability. Figure 1.2 is a 
schematic chart of the nuclides adapted from Norman and Schramm 
(1979) illustrating the neutron capture path of the s-process. The hor- 
izontal axis is the neutron number, N, and the vertical axis is the atomic 
number, Z. Stable nuclei are indicated by dots. Figure 1.3 is a detail of 
the chart in the Sn-Ce region where the stable nuclides are now indi- 
cated by a box with the mass number, A, inside. If we start at na Sn, the 
s-process path proceeds by neutron capture through the heavier iso- 
topes of Sn until reaching lzl Sn which, unstable with a half life of 2.8 
hours, undergoes beta decay and becomes 121 Sb before capturing 
another neutron. As a result, 122 Sn and I24 Sn are not formed in the clas- 
sic s-process. The path continues in this same fashion along the valley 
of beta stability, bypassing the isolated neutron rich isotopes such as 
3ZS Te and l34 Xe. 


The classic r-process is at the other extreme. Here the flux of neu- 
trons is extremely large with the neutron capture times assumed Lo be 
very much less than the beta decay times of the nuclides involved. As a 
result, the path of the r-process progresses through extremely neutron 
rich nuclei far from the valley of beta stability. The location of the path 
is determined by (n, 7 )~ ( 7 ,n) reactions which prevent further neutron 
addition when the binding energy of the extra neutron is insufficient to 
prevent dissociation by the ambient photon gas (the neutron drip Line). 
In a more generalized r-process (Blake and Schramm 1976), the "path'' is 
determined by the place where the beta decay times of the increasing iy 
neutron rich nuclei become comparable with the neutron capture times. 


Figure 1.2 


A schematic chart of the nuclides adapted from Norman and 
Schramm (1979) showing the paths of the r-process and the s- 
process. Neutron number is along Ur- horizontal axis and 
atomic number is along the vertical. Stable nuclei are indicated 
by dots. Also shown are the neutron and proton magic numbers. 
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Figure 1.3 


A detail of the chart of the nuclides in the soSn - M Ce region. 
The vertical axis is the atomic number (element) and the hor- 
izontal axis is the neutron number. Mass numbers are indicated 
for the stable nuclides. The light diagonal lines are meant to 
suggest the beta decay of the results of r-process nucleosyn- 
thesis. 
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although use of the word path here does not imply tha^ only one unique 
capture sequence is allowed. One possible r-proeess path is indicated in 
Figure 1.2 and typically involves nuclides some 10 or more neutrons 
richer than the valley of beta stability. It is only after the neutron flux 
ceases that the nuclei beta decay back to the valley of beta stability. In 
Figure 1.3 this is indicated by the diagonal lines terminating at the first 
stable nuclide encountered. Note that, as a result of this, nuclei such as 
IZ4 Te and 130 Xe cannot be produced by the r-process because they are 
shielded from the beta decay of the r-process path by l34 Sn and 130 Te, 
respectively. Such nuclei are culled s-only nuclei. 

The proton rich nuclei to the left of the valley of beta stability can- 
not be produced by either the r- or the s-process. Another process, the 
p-proeess, a proton capture or neutron removal process, is invoked to 
explain their abundances (Arnett 1973 discusses the p-process briefly in 
reference to explosive nucleosynthesis). However, since these isotopes 
are much less abundant (generally 1-2 orders of magnitude) than s- and 
reprocess isotopes, we will ignore contributions due to the p-process in 
what follows. 

The determination of the abundances resulting from the s-process 
involves knowledge of the neutron capture cross sections of all nuclei 
along the s-process path. In the s-process, cr A N A , where N A is the s- 
proeess abundance of a nuclide along the path with mass number A and 
neutron capture cross section cr A , is a smoothly varying function of A, 
and, at least locally, the relative abundances of nuclides are inversely 
proportional to their neutron capture cross sections. Nuclei having a 
magic number of neutrons (and to a lesser extent a magic number of 
protons) have small neutron capture cross sections and therefore large 
s-process abundances. This is the ease for ,38 Ba which has a magic 
number of neutrons and also for 50 Sn which has a magic number of pro- 
tons. The validity of the s-process model can be tested by examining its 


ability to reproduce the known abundances of the s-only isotopes. 
Results of Kappeler et al. (1982), for example, show that it is able to do 
this quite well with the correct choice of neutron exposure values, within 
the uncertainties in the measured cross sections and abundances. 

The situation for the r-process is not as well established. No site 
responsible for r-process synthesis has as yet been agreed upon, 
although attempts have been made to place limits on the temperature 
and density of the syn thesis region (e.g. Norman and Schramm 1979) 
both of which affect the location of the r-procass path. The problem is 
further complicated by the need for neutron capture and beta decay 
rates for nuclei far from the valley of beta stability, the large majority of 
which have never been synthesized in the laboratory. The necessary 
information is obtained by extrapolation using various nuclear models 
with order of magnitude differences between beta decay rates from 
different models. Nevertheless, the r-process is similar to the s-process 
in that the production of neutron magic number nuclei along its path is 
favored because their longer decay times result in a buildup at these 
magic numbers. However, because the r-process path is displaced from 
the valley of beta stability during the neutron exposure, the subsequent 
beta decay results in a shifting of the r-process pciaks to lower atomic 
numbers than those of the s-process, as can be deduced from an exami- 
nation of Figure 1.2. The actual amount of the shift depends on the 
location assumed for the r-process which in turn is a function of neu- 
tron density, ambient temperature, and beta decay rates (see e.g. Nor- 
man and Schramm 1979, Schramm 1973, and Cowan et al. 1983). 

The results of r-process calculations can be compared with the "r- 
only" isotopes (although the existence of r-only isotopes may be only an 
approximation if the s-proeess occurs in flashes, see Cameron 1902a). 
What is usually done however, is to decompose a given abundance distri- 
bution into r- and s-process contributions by subtracting the s-process 


theory results from the full distribution. This is done by matching the 
s-process model to the s-only isotopes and then using the a k N A curve to 
deterrriine the s-process contribution to the mixed isotopes. Figure 1.4 
shows a decomposition of the solar system abundances in the Sn-Ce 
region into r and s process contributions. As explained in more detail in 
Chapter 4, this was derived using the s-process of Kappeler et al. (1982) 
and subtracting it from the solar system abundances of Anders and 
Ebihara (1982) on an isotope by isotope basis. Other decompositions, 
such as Cameron 1982a, differ in details but the general structure 
remains similar. The main features to note here are that, for the s~ 
process contribution, 5 oSn and 56 la, the magic number nuclei, are 
significantly more abundant than 52 Te and ^Xe. For the r-process, the 
situation is reversed with Te and Xe being the dominant elements. Addi- 
tionally, 5fl Ce can be seen to be primarily an s-process element. These 
rather gross differences make this region a sensitive indicator of the 
relative amount of r- and s-process material present. 

Early measurements of the ultra-heavy cosmic rays in a higher 
charge range indicated the presence of an r-process contribution. 
These experiments were primarily detectors of the passive type consist- 
ing of large areas of plastic track detector and/or photographic emi'.!- 
sion flown on balloons and later recovered. Charge identification was 
made by examining the damage trails left by the particle where it 
penetrated the material. (The ultra-heavy component of the cosmic rays 
was first detected by means of "fossil" tracks in meteorites (Fleischer et 
al. 1967)). The need for hand scanning of the detector material resulted 
in attention being focused on the high charge, rarer nuclei for which the 
detection efficiencies were near unity. The s-proeess is unable to pro- 
duce elements heavier than because of the lack of stable nuclei 

between 83 @i and aoTh. Therefore, detection of any nuclei with charge 
greater than 03 Bi, in particular the actinides 9Q Th and B 2 U, would be 
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Figurs 1.4 

A decomposition of the solar system material into r- and s- 
process contributions using the solar system abundances of 
Anders and Ebihara (1982) and the s-process of Kappeler et a 1. 
(1982) (See Chapter 4). The solar system abundance of u Si is 
defined to be 10®. 


ABUNDANCE (Si ■ 10 ) 





evidence of r-process contributions. Several such particles were 
reported both in early balloon results (Fowler 1977 gives a brief review) 
and in measurements made from long term exposures of Lexan on board 
Skylab (Shirk and Price 1978). However, the charge resolution and/or 
exposures were such that only element groups, and not individual 
charges, could be studied. 

Electronic detectors, with their ability to process greater amounts 
of data and their promising charge resolution, were used to study the 
region just above iron (Israel et al. 1979). The results agreed better with 
a solar system mixture of r- and s-proeess material than with either 
component separately. However, because of the combination of small 
instrument geometry and short exposure times on balloons, statistics 
were insufficient to measure abundances above charge 40. Only with the 
advent of the long exposure times available on satellites could the rarer 
elements, such as Sn and Ba, be studied. Results have recently become 
available from two electronic satellite experiments to detect the ultra- 
heavy cosmic rays: the University of Bristol experiment on the Ariel VI 
satellite and the Heavy Nuclei Experiment on the HEAO-3 satellite. 

This thesis will report on measurements of the abundances of the 
even charge elements in the region from soSn to 5 aCe based on a first 
study of the data from 440 days of operation of the HEAO-3 Heavy Nuclei 
Experiment. These results are inconsistent with a pure r-process source 
hut do not rule out the presence a solar system type mixture and are 
therefore consistent with measurements made by the same instrument 
in the charge range 26 :£ Z s 40 (Binns et at. 198 lb) and in the actinide 
region (Binns et al 1982). 
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II. THE INSTRUMENT 


The Heavy Nuclei Experiment on the third High Energy Astronomy 
Observatory (HEAO-3) satellite was designed to measure the elemental 
composition of the ultra-heavy component of the cosmic radiation. The 
experiment was the result of a collaboration of three research institu- 
tions: the California Institute of Technology, Washington University at St. 
Louis, and the University of Minnesota at Minneapolis. A description of 
the instrument, its electronics, and its design considerations can be 
found in Binns e£ of. (1961a). Here I will give a summary description of 
the instrument emphasizing those features which are relevant to the 
present results. 


Because the flux of the ultra-heavy cosmic radiation is many orders 
of magnitude below that of the more abundant nuclei of carbon, oxygen, 
and iron, any instrument designed to measure its elemental composition 
w: h reasonable statistical accuracy must, of necessity, have a large col- 
lecting power (geometry factor) and/or long exposure times. As a 
result, one of the most striking characteristics of the Heavy Nuclei 
Experiment instrument is its size. Whereas most cosmic ray detectors 
designed for use on spacecraft to measure abundances in the lower 
charge ranges have areas in the neighborhood of 500 mm 2 and geometry 
factors of less than or equal to 1 cm 2 sr, this experiment has an area of 
about 2 m 2 and a toteil geometry factor of 59,000 cm? sr. 

The determination of the charge of a particle incident on the detec- 
tor is by means af the d'E/dx- Cerenkov technique. In this method the 
particle passes through both an ionization chamber, which measures the 
amount of energy deposited in the detector through interactions with 
atomic electrons, and a Cerenkov counter, which measures the amount 
of Cerenkov light emitted as the particle penetrates the radiator at a 
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velocity greater than that of light in the radiator material. The response 
of both of these detectors is a function of the speed of the particle 
(divided by the speed of light), p, and its charge, Z. The ionization signal, 
I, can be •written approximately as (Janni 1866) 


zi 


— p 2 + k 


( 2 . 1 ) 


while for the Cerenkov signal, C, 


C pcZ 3 



( 2 . 2 ) 


where 7 = (1 — k is approximately a constant, and n is the index of 

refraction of the Cerenkov radiator material. Since the cosmic rays 
have a spectrum of incident energies and charges, both of these meas- 
urements are, in the absence of other information, necessary for the 
determination of the particle's charge. 

Figure 2.1 is plot of the square root of the ionization signal per unit 
pathlength, ZI, versus the square root of the Cerenkov signal per unit 
patMength, ZC, Shown on this plot is a family of curves representing the 
response of the HEAO instrument to the elements wit ft even charge in 
the iron region. Each curve represents the response to a different 
charge with distance along the curve parametrized by the particle 
energy. The minimum in ZI for each curve is at approximately 2 
GeV/ nucleon with the increase in ZI at high energies, the relativistic 
rise, the result of the logarithmic term in (2.1). For both the Cerenkov 
and the ionization response we have used the the Z-sqnared scaling 
shown above and in the lower energy regime slowing down of the particle 
within the detector has been taken into account (Krombel 1960). One 
should note that for high enough energies the Cerenkov response alone 
is sufficient to determine the charge. 
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Figure 2.1 


A plot of the square root of the ion chamber signal per unit 
pathlength, ZI, versus the square root of the Cerenkov chamber 
signal per unit pathlength, ZC. The curves on the plot show the 
expected response of the Heavy Nuclei Experiment to the even 
charge elements from uSi to ^Se. The normalization of ZC is 
such that a charge 26 (iron) particle will result in a ZC of 
26. ZI is normalized so that its minimum value for iron is 26. 
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Figure 2.2 is a cross section view of the experiment. The instrument 
consists of six parallel plate ionization chambers (IC), four pairs of mul- 
tiwire ionization hodoscopes (H0D0), and a dual radiator Cerenkov 
detector (RADI and RAD2). The arrangement is symmetric about the 
detector midplane (the X-Z plane) with three ion chambers, two sets of 
hodoscope planes, and a Cerenkov radiator on either side. Each "half" 
of the instrument consists of a sealed pressure vessel containing the 
three ion chambers and two hodoscopes. The Cerenkov radiator is 
mounted on the outside of the pressure vessel and the module so formed 
is attached to its mate with the Cerenkov radiators facing each other. 
The radiators do not touch but are separated by a vacuum and the 
Cerenkov light box is formed by enclosing this volume in a light-tight 
seal. Both radiators are viewed by the same set of eight photo multiplier 
tubes. The instrument and each of the modules has the shape of a rec- 
tangular parallelepiped which was dictated by the satellite geometry. 

Each pressure vessel is a sealed unit in. order to obviate the neces- 
sity for an external onboard gas supply. P-10 gas has been used, which 
is a mixture of 90% argon and 10% methane with a trace (approx 0.5%) of 
helium added for leak detection. The gas pressure is a nominal 838 torr 
at 20° C. The instrument windows consist of aluminum honeycomb with 
a thickness of 8.9 cm and a mean areal density of 1.2 g cm -2 . Aluminum 
honeycomb was chosen for the window material because it combines the 
features of high strength and low density thereby minimizing the 
amount of fragmentation which occurs as a particle penetrates the 
instrument. The ultra-heavy cosmic rays are particularly sensitive to 
this because of their small fragmentation pathlengths. 

The ionization chambers are of dual-gap design with the anode 
mounted midway between the two cathodes. The cathodes are shared by 
adjacent ion chambers and/or the adjacent hodoscope planes. The elec- 
trodes are made of aluminum screenwire and the ion chambers have an 


Figure 2.2 


A cross sectional view of the Heavy Nuclei Experiment on HEAO- 
3 showing the ion chambers (1C), the Cerenkov radiator (RAD 1 
and RAD 2), and the hodoseope planes (HODO). The aluminum 
honeycomb windows are suggested by 4 he dashed lines showing 
the division of the instrument into two separate pressure 
vessels. 












anode to cathode spacing of 6.9 cm and are operated at -1000 V. This 
choice of operating voltage was made to minimize both electron collec- 
tion losses and variation of electron drift velocity with small variations 
in pressure or voltage as detailed in Binns at al. (1901a), Each ion 
chamber is separately pulse height analyzed and all six pulse heights 
recorded for each accepted particle. 

The trajectory of the incident particle is determined by the use of 
discrete wire ionization hodoscopes. Each hodcjcope layer consists of a 
an anode composed of 0.025 cm diameter parallel stainless steel wires 
with a center to center spacing of 1 cm The anode layer is midway 
between two screen wire cathodes with an anode to cathode spacing 
again of 1 cm. The operating voltage is -1000 V and the diameter of the 
wires is such that the hodoscopes are operated in the ionization mode 
and no gas amplification takes place. Each anode wire has its own 
charge sensitive preamplifier and discriminator. For each event, the 
instrument records the discriminator state of up to 16 wires in each of 
the 8 hodoscope layers in the form of two address /patterns. The first 
address /pattern consists of the address of the first (lowest address) wire 
fired and the discriminator state of the subsequent 7 wires. The second 
pattern has the same format as the first and is used when more than 8 
wires have fired. 

The Cerenkov counter portion of the detector consists of two sheets 
of 0.47 cm thick Pilot 425 viewed by eight photomultiplier tubes. Each 
sheet has been sandblasted in order to improve uniformity of response 
and the face against the pressure vessel has been painted with white 
paint as has the interior of the light box. The two radiators are 
separated by a distance of 24.7 cm. The eight photomultiplier tubes are 
arranged in pairs with one pair at each corner of the rectangle forming 
the light box. Each photomultiplier tube is separately pulse height 
analyzed and the eight values are recorded for each event. The index of 


refraction of the Pilot 425 for the wavelengths of interest is approxi- 
mately 1.52 which corresponds to a threshold energy of approximately 
300 MeV/amu. 


The criteria for acceptance of an event for analysis, in the normal 
mode of operation, are triggering of one wire or more in both the X and 
Z planes in at least two of the four hodoscopes and the firing of at least 
two of the seven charge measuring detectors, i.e. the six ion chambers 
and the Cerenkov detector. These particular requirements were meant 
to ensure that one is able to make a trajectory determination for the 
event but do not ensure that the event has both a Cerenkov and an ion 


signal. A perfectly valid event could conceivably pass through two hodo- 


scopes and miss both Cerenkov radiators (e.g. triggering the Xl-Zl HODO 
and the X2-Z2 HODO in Figure 2.2). One can still assign a charge to this 


kind of event, albeit with reduced confidence, using geomagnetic cutoff 
data to restrict the possible energy rarge. However all the events used 
in this analysis were required to have a Cerenkov signal. 


The minimum charge particle needed to trigger the hodc^ccpe 
depends on the the angle of incidence of the particle and its position 
with respect to the anode wire. For the typical trajectory the path con- 
tributing to one anode wire is about 2 cm and the minimum charge par- 
ticle that will fire the discriminator on that wire is about 11.5 charge 
units (all charges are in units of the proton charge). A vertical trajectory 
midway between two anodes, the worst ease, has a path oi 1 cm contri- 
buting to each anode and requires a particle with charge of approxi- 
mately 16.2 to trigger either discriminator. The Cerenkov trigger, C2, 
consists of the "or" of two photomultiplier tubes not in the same corner 
of the instrument as shown in Figure 2.3. The discriminator for each 
photomultiplier tube is set at 16% of a perpendicular £ = 1 charge 25 
(iron) nucleus or the full signal of a perpendicular, Z ~ 10, /S ~ 1 nucleus. 
The actual charge threshold is e function of incident angle because 


- 25 - 


Figure 2.3 

A schematic drawing of the Cerenkov radiator layout showing 
placement and labeling of individual photomultiplier tubes along 
with the Cerenkov triggering criterion. 



particles with a larger angle to the normal have, in general, a greater 
pathlength in the detector and hence a larger signal. The ion chambers 
each have three different discriminator levels: LLP, ILD, and HLD. These 
three levels are jet to correspond to vertically incident charges 16.1, 
33.3, and 43.7 respectively. The LLD is the trigger used lor event 
analysis while the other discriminators are used for determining event 
''priority". The presence of a relatively large flux of iron nuclei results in 
a significant amount of dead time for non-priority events due to the data 
transmission rate of 128 events every 40.96 seconds. The priority sys- 
tem was instituted to ensure that the high charge events are recorded 
with essentially 100% efficiency by allowing a priority event to write over 
any non-priority event being held for transmission. An event is desig- 
nated priority if any HLD fires or if any ILD and C2 fire. These ensure 
that any event with charge greater them 44 is recorded as well as any 
event with charge greater than 33 and an energy greater than about 350 
MeV/ nucleon. 

Figure 2.4 is a schematic view of the satellite indicating the place- 
ment and configuration of the Heavy Nuclei Experiment. The instrument 
windows sure exposed on either side of the spacecraft vith the instru- 
ment itself oriented along the spacecraft's Y axis. The rather loose 
event acceptance criteria allow for the analysis of particles which did 
not enter through the windows. Such "sidewall" events have had to 
penetrate larger amounts of material than the window events, typically 
2-3 g cm -2 , due to the satellite body and electronic packages mounted to 
the outside of the instrument. However, allowing such events to be 
analyzed results in approximately a threefold increase in our geometry 
factor over using window events exclusively. 

The spacecraft was launched into a circular orbit with initial alti- 
tude of 496 km and an inclination of 43, .6° on 1979 September 20. The 
usual mode of operation was for the spacecraft Z axis (the solar par el 


Figure 2.4 


A schematic view of the HEAO-3 satellite showing the placement 
of the Heavy Nuclei Experiment (Linder 1979). The spacecraft is 
normally oriented with its Z axis (and solar panels) pointed 
towards the sun. It spins about the Z axis with a 20 minute 
period. 
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axis) to be pointed toward the sun with a rotation period of 20 min 
around this axis. As a result, the instrument does not have a fixed 
orientation with respect to the earth. To accommodate the other exper- 
iments on board, there were several periods when the Z axis was 
directed towards a point some 30° away from the sun ("offset scan" 
mode). Star sensors on the satellite allow for post-facto attitude deter- 
mination to better than 0.5°. 


The altitude and inclination of the orbit are such that the space- 
craft will pass though a region of the geomagnetic field known as the 
South Atlantic Anomaly (SAA) during many of its orbits. The 
configuration of the geomagnetic field is such that a large flux of 
trapped protons are encountered here with the flux being high enough 
that the instrument response is severly degraded by accidental coin- 
cidences, possible baseline shifts, and other effects. As a result data 
taken during an SAA passage are extremely unreliable without special 
processing and are not used in this study. 

The Instrument operated until 1981 May 29 with the latter part o; 
the mission characterized by a degradation in response in several of the 
ion chambers (for that time period, the ion chambers showing the degra- 
dation were not used in the present analysis). We report here on the 
results from approximately 440 days of operation. Appendix A lists the 
time periods used. 
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H. THE DATA 


3.1. O v e r view 

Because of the large variation in the flux of the ultra-heavy cosmic 
rays (a decrease of at least five orders of magnitude from iron to 
uranium) and the desire to establish a reference point in a previously 
investigated charge range, a large number of iron, and other lower 
charge nuclei, are present in the data telemetered back to earth from 
the Heavy Nuclei Experiment on the HEAO-3 satellite. Typically some 
180,000 events per day were recorded by the instrument. The presence 
of such a large flux of particles allows for in-flight calibrations and map- 
ping of instrument response. However, the sheer volume of data creates 
problems in data processing and handling, especially when we consider 
the amount of sorting that must be done to select out the less than 150 
events in the charge range of interest here. 

In order to facilitate analysis the data processing was divided into 
several stages, each more selective than the preceding. The first stage, 
or pass through the data, converts all events from their raw encoded 
form into one which is more easily interpreted. In the second pass, a 
high charge subset is selected from the output of the first pass and the 
charge estimates are further refined by more sophisticated processing. 
Events from this greatly reduced subset were classified on the basis of 
their probable energy and separated into categories of different 
expected charge resolution. The various categories were then examined 
using the technique of modulo 2 superposition of even element peaks to 
select those actually having the best resolution. Those selected were 
added together to yield the final results presented in Section 3.6. In this 
chapter, we will describe in more detail the processing steps outlined 
above. 


3.2. First Pass Analysis 


In the first pass at data analysis no attempt was made to be selec- 
tive but instead as many events as possible were converted from their 
raw form, as telemetered from the satellite, into more physically mean- 
ingful quantities. This processing was done by the library generator pro- 
gram, LIBGEN (Garrard 1979a) which converted production (raw data) 
tapes, as received from Goddard Space Flight Center, into library tapes. 
The data from the instrument fill one 2400 foot, 1600 BPI production 
tape per day of operation and each Library tape contains the analyzed 
data from one production tape. 

The LIBGEN program is a rather extensive algorithm which does the 
following: 

1. Converts the wire patterns from the instrument hodoseopes into 
trajectories using known dimensions for the wire and hodoscope 
spacings. This step results in the first major classification of parti- 
cle types depending on the accuracy with which a trajectory can be 
constructed. A list is maintained of wires in each layer for which the 
firing rates have been determined to be excessively high or low 
based on examinations of daily rate data. This information is used 
to eliminate "missing" or "extra" wires from the hodoscope patterns. 
For those events which have patterns with no "missing" or "extra" 
wires in 2 or more hodoscope planes (in both the X and Z coordi- 
nates of that plane), the center of each pattern is used in a least 
squares fit to a straight line to obtain the particle trajectory. Those 
events having inconsistencies in the wire patterns which preclude 
construction of reliable trajectories are classified separately from 
the good trajectory events for future consideration. 


Flags those events which have data quality or logic problems, i.e. 
parity or telemetry errors, nominally repeated events which do not 
match their predecessor, etc. Also flagged are those events for 
which trajectories were able to be constructed, but which still had 
trajectory inconsistencies such as more wires triggered than is con- 
sistent with the calculated trajectory, an ion chamber which is cn 
the trajectory but which does not have a signal, a large x 2 on the 
fitted trajectory, etc. 

Converts raw ion chamber pulse heights, from the pseudo- 
logarithmic pulse height analyzers (PHA), to signal in femto- 
coulorribs. 

Converts the raw photomultiplier tube pulse heights, again from a 
pseudo-logarithmic PHA, into signal in volts. 

Normalizes both ion and Cerenkov signals to signal per unit path- 
length using the trajectory information and measurements made 
prior to launch of the thickness of the Cerenkov radiators. Each 
radiator has its own thickness map which is used in the calculation. 


Extracts and processes the satellite attitude and orbit data for each 
event which are necessary to relate the particle trajectory in the 
spacecraft frame to an external frame of reference. 

Makes initial estimates of the particle’s charge using Cerenkov and 
ion chamber data. These estimates are based on the fact that the 
ratio of the Cerenkov signal to the ion signal is, to first order, 
independent of particle charge with both being proportional to Z 2 in 
this approximation. Thus, referring to Figure 2.1, one can derive a 
charge independent function of ZC/ZI, using iron data, which will 
correct ZI by an amount determined by ZC/ZI and allow an energy 
independent estimate of the charge to be made (Israel 1900 and 
private communication). Such a function can be applied to the 
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average ion chamber signal from either module. If either of these 
estimates is greater than 30.5, the event is flagged as a "high 
charge” event. For 54% of the time period reported on here, ZC 
alone was also used to select high charge events. Because of subse- 
quent ion chamber and Cerenkov agreement criteria and the dis- 
tance of the threshold from the present charge region, this distinc- 
tion is unimportant for this analysis. 


3.3. High Charge Subset 

The first pass does not decrease the volume of the data but does do 
a preliminary classification and processing of the events. In this next 
pass only those events which had some indication of having a charge 
greater than 30.5, i.e. the "high charge” events mentioned above, were 
retained. In addition, more refined processing was done to obtain better 
charge estimates. 

Because of the large surface area of the detector, it was necessary 
to make corrections for nonuniformities in detector response over the 
active area of the instrument. For this purpose the large flux of iron 
nuclei was useful in deriving a two dimensional response map for each 
ion chamber and for the Cerenkov detector. In mapping the ion 
chambers, variations in the ionization signal caused by differences in 
particle energies were reduced by choosing a subset of the iron data in 
the minimum ionizing region. All the particles in this subset should have 
the same ionization signal and were binned according to their position 
within the chamber under study. The resulting maps are uniform in the 
central area (more than 44 cm from any wall) with variations of less 
than 0.1%. The response falls off linearly to 0.98 of the central value as 
one approaches to within 8 cm of a wall. For events used in this 
analysis, we ignore all chambers for which the particle trajectory indi- 
cates passage closer than 8 cm to a wall within that chamber. These 


maps were constructed at Washington University and each ion chamber 
has its own position dependent correction based on its response map. 
The in-flight iron data were also used to normalize the chambers to one 
another. 


The Cerenkov chamber map, on the other hand, has a much larger 
variation over its surface. Response maps were made for each individual 
photomultiplier tube and for the mean of all eight tubes. Only the eight 
tube map was used in the response corrections, however, with the single 
tube maps being used solely for the determination of consistency cri- 
teria (see below). The largest light collection non- uniformities are the 
result of "hot spots" in front of each PMT. In the central region of the 
eight tube map, more than 25 cm from any PMT, the typical gradients 
are less than 0.2%/cm with variations of about 10% over the entire cen- 
tral region. The corrections in front of a PMT however can be as large as 
0.72 with gradients of 1.5 %/ciel The iron data were again used to con- 
struct the Cerenkov maps. This work was done at the University of Min- 
nesota. 

In addition to the mapping corrections, there was a correction 
applied for time variations in response. Observation of the iron peak 
position showed that the ion chamber time variation was less than 0.2%. 
The individual PMTs however showed a variation which correlated well 
with temperature and was of the order of 0.7 %/°C. However, the gain 
variations were slow enough, usually well under 0.5 %/daj r , that sufficient 
numbers of iron particles were collected to allow corrections to be made 
on a daily basis, to each PMT, to 0.15 %. Both the ion chamber and the 
photomultiplier tube response were monitored on a daily basis for the 
duration of the flight. In addition, corrections were made for differences 
in the individual PMT gains so that all tubes were weighted equally in the 
final charge determination. 


Beth the mapping and the corrections for daily variations in the 
PMTs were done in this second pass. Additionally, a more sophisticated 
algorithm was applied to the particles for which the LIBGEN program was 
unable to determine a trajectory in order to assign, by deletion of bad 
hodoscope planes etc., a trajectory and a charge to these events. How- 
ever, neither these events nor those flagged as having residual trajec- 
tory inconsistencies were used in the data set presented here. Elim- 
inated before the next analysis stage were those events previously 
flagged as having data quality or logic problems. In addition, preliminary 
consistency checks were made using the multiparameter nature of the 
instrument in order to remove obviously anomalous events before 
further processing. 

The first of these selections was based on agreement between the 
measurements made by the individual photomultiplier tubes for a given 
event. Figure 3.1 exhibits histograms, for the high charge data set, of 
the ratio of single tube to mean of all eight both for tubes which are in 
the same quadrant as the panicle and for tubes which are in the oppo- 
site X-half of the Cerenkov box. No single tube mapping has been done 
for this selection and the width of the distributions is dominated by the 
mapping variations. The difference in the peak location between the two 
histograms indicates that the near quadrant tubes have a signal that is 
some 1 .3 times the average while the distant tubes have only 75% of the 
average signal. By using the individual tube maps, limits on the allow- 
able deviation of a single tube from the mean of eight tubes was deter- 
mined. The agreement criteria were based on the location of the tube 
with respect to the particle’s position at the midplane of the Cerenkov 
box. Only a gross determination of location was used: whether the tube 
was in the same quadrant (or half) of the radiator as the particle. This 
selection helps to eliminate errors in the Cerenkov charge due to single 
tube errors. By far the most common tube error is the presence of one 
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Figure 3.1 

Histograms of the ratio of individual Cerenkov photomultiplier 
tubes to the mean of all eight tubes. The eight tube average has 
been corrected for position variations in response. The indivi- 
dual tube values have not. The upper histogram is for tubes 
determined to be in the same quadrant as the particle based on 
its position at the center of the Cerenkov box (the "closest" 
tubes). The lower histogram is the same for tubes in the oppo- 
site X-half of the Cerenkov box (the "farthest" tubes). Also indi- 
cated are the agreement limits used in the consistency selec- 
tions. "Single tube hits" have been eliminated from these plots. 





tube which dominates the sum ef eight. That these events are the result 
of the particle actually hitting the photomultiplier tube in question is 
borne out by the particle trajectory which almost always indicates that 
it passed through the Cerenkov box in a corner. These spurious "high 
charge” events constitute some 00% of the high charge data set. The 
"single tube hits" have been removed from Figure 3.1 by eliminating 
those events having a ratio greater than 6.0. (For the case where one 
tube dominates the average, we would expect its ratio to the average to 
be approximately 0.) The location of the consistency cuts are indicated 
on the histograms. The accepted values were between 0.6 and 6.0 for the 
same quadrant tubes, 0.5 and 2.5 for same X-half but not same qua- 
drant, and 0.4 and 1.6 for opposite X-half tubes (Israel 1979). 

A second selection was made on the agreement between ion 
chambers within a module. The amount of matter between the chambers 
in one module is small enough (<0.1 g cm -2 ) that there should be no 
appreciable energy loss from one chamber to the next. Figure 3.2 shows 
a histogram of the dispersion between chambers within a single module 
(the module containing ion chambers 1, 2, and 3 in Figure 2.2). The his- 
togram includes both two and three chamber events. The quantity plot- 
ted is 


n 2 I ^ kyg I 

FDZ=l-^~ (3.1) 

* i'est 

where I* is the signal from ion chamber i, I ttVg is the mean ion signal 
within that module, Z e#t is the charge estimate (of the type discussed in 
section 3.2) based on the module, and the sum. over i (and I OV g) includes 
only the n ion chambers in the module with a valid signal. This value, in 
the absence of correlated errors between chambers, corresponds to the 
uncertainty in the charge estimate expected from using the ion 
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Flgure 3.2 

Histogram showing the degree of agreement between ion 
chambers within a single module. The horizontal axis is a meas- 
ure of the dispersion between chambers, in charge units (see 
text). Also shown is the consistency selection used. 



-42 - 


chambers since, because of the Z-squared dependence of the signal, I, 
I * Z 2 and thus 2 ~ = — The amou n t of dispersion indicated by the 

mode of this plot however, approximately 0.3 charge units (measured in 
terms of the proton charge), is lower than would be expected on the 
basis of the observed ion chamber charge resolution, and is probably the 
result of correlations between chambers within a module. In order to 
eliminate only those events having obvious ion chamber inconsistencies, 
we required FDZ < 2.25 change units. 

As mentioned in the previous chapter, events which occurred during 
passage through the South Atlantic Anomaly (SAA) were also excluded 
from consideration. The location of the SAA was determined both on the 
basis of orbital position and the "singles" rates in the ion chambers, i.e. 
the number of firings of the ion chambers as indicated by the number of 
signals above the LLD, whether or not the event met the coincidence 
requirement. 


3.4. Energy Selections 


For this first study, because it has better resolution, we have used 
the Cerenkov chamber alone as the determinant of charge with the ion 
chamber information used only for consistency checks and in the 
classification scheme. As a result, it is necessary to accept only those 
particles having energy great enough that the Cerenkov response has 
reached its plateau, This can be done by using the earth’s magnetic 


field as an rigidity filter (Rigidity, R, being defined as 


SR. 

Ze 


where p is the 


particle momentum, e is the speed of light, and Z is the charge of the 
particle in units of the proton charge, e). In the earth’s field, particles 
arriving from any given direction must have a rigidity greater than some 
critical rigidity which is a function of the direction of arrival and posi- 
tion of observation in the field. The theory of motion in the earth's 
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ra&gnetie field was first developed by Stdrmer (e.g., Stdrmer, 1955) by 
modeling the geomagnetic field as a dipole. The axial symmetry of the 
dipole results in a constant of the motion which allows derivation of the 


following equation for the Stdrmer cutoff rigidity, R 0 , as a function of 


observation position in the dipole and direction of viewing. 


p _ M. 


1 — [l — cosy cos s a| 


cosy eosA 


(3.2) 


where M is the magnitude of the dipole moment, r is the radius from the 
center of the dipole, y is the angle that the particle’s trajectory makes 
with the west (a particle traveling due west has a y of zero) and A is the 
latitude of the observation point. Although the geomagnetic field is not 
truly a dipole and the determination of the exact Re involves compli- 
cated calculations for tracing the particle's assumed trajectory back- 
wards, we can still define a local magnetic west and local magnetic lati- 
tude using the value of the magnetic B field and the Me Twain L value at 
the point in question. Then using the dipole approximations (Roederer 
1970) that 


B 2 = ^—[4 — 3cos 2 A ] 

and 


L - 



(3.3) 


(3.4) 


(with r e the radius of the earth) we can eliminate A and estimate the 
cutoff rigidity for the direction of arrival of the particle (Garrard 1979b). 
The values of B and L used were calculated at Goddard Space Might 
Center from standard reference geomagnetic fields ’vhieh approximate 
the earth’s actual field by a spherical harmonic expansion. Figure 3.3 is 


a plot of ZI versus ZC for data in the sub-iron region Loth with and 
without a rigidity selection applied. It can be seen that, although the 
theory used is an approximation, the rigidity cut does succeed in exclud- 
ing the majority of the low energy particles. A cutoff rigidity of 8 GV was 
used for this analysis. For a0 Sn (typical atomic mass 120 amu), this 
corresponds to an energy of 2.53 GeV/nuelecii and a Cerenkov signal 
which is 94% of the (i = 1 value. 

The determination of a unique R*. involves knowledge not only of the 
particle trajectory, but also of the direction of motion along the trajec- 
tory. Changing y by 180° changes R,,. Since the Heavy Nuclei Experi- 
ment does not have a device for determining direction of motion along 
the trajectory, two possible R c are assigned to each eve/ 1, corresponding 
to the two possible directions. For some of these particles, only one of 
these two directions is permissible because, if traced Backwards in the 
other direction, the trajectory intersects the earth. Determination of 
these earth shadowing directions was done empirically by using data 
provided by the Danish-French cosmic ray isotope experiment also on 
the HEAO-3 satellite (Lund and Westergaard, private communication; 
Garrard and Ennis 1980a 8c b). Having a time of flight device, they were 
able to map out forbidden directions in the sky. In practice, use of the 
forbidden directions is necessary for only 12% of the high rigidity data 
selected for use in the Sn-Ce region. The remainder of the events have 
both Rc greater than 8 GV. 

Since the high cutoff particles constitute only 40% of the data, in 
order to increase the number of events in the sample we have also 
analyzed particles which were chosen to be high energy on the basis of 
their ratio of ZC to ZI. Figure 3.4 is a crossplot of ZC/ZI versus ZC for 
two days of selected iron data which illustrates the rationale behind this 
selection. The curves in Figure 3.4 are the same as the element tracks 
seen in Figures 3.3 and 2.2 and, again, distance along the curve can be 


Figure 3.3 


Crossplots of the ion chamber charge estimate versus the 
Cerenkov charge estimate for two days of data in the iron and 
sub-iron charge region. Only events with two estimates of the 
charge based on the ion chambers, i.e. two module events, which 
agree to within approximately 7.5% were used. The upper plot 
has only particles with rigidities above 8 GV, chosen using the 
selection in the text. The lower has no rigidity selection. (The 
lack of events at ZC<7 is an artifact of the selection program.) 
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Figure 3.4 

Crossplot of the ratio of ZC to ZI as a function of ZC for two days 
of iron data. This plot illustrates how the selection on ZC/ZI 
(see text) chooses only the "tip*' of the response curve enabling 
the use of ZC as a valid charge estimate. Note that silicon 
(Z=14) may be affected by threshold effects. 





parametrized by the particle energy. (In Figures 3.3 and 2.2, curves of 
constant ZC/Zf are straight lines through the origin with slope equal to 
the ratio.) Events on Figure 3.4 with the highest values of ZC/ZI, at the 
tip of the curve, correspond to energies in the neighborhood of 5 
GeV/nucleon with lower energies lying to the lower left and higher ener- 
gies to the lower right. By selecting on ZC/ZI, we can eliminate the low 
energy particles of each element which would otherwise contaminate the 
charge peaks below it due to the energy dependence of the Cerenkov sig- 
nal. The value used in the selection, 0.964, was chosen by varying the 
ZC/ZI requirement on the iron data in order to obtain the best combina- 
tion of statistics and resolution. 


3.5. Event Selections 

Having selected events for which the Cerenkov charge estimate, ZC, 
should be a valid measure of the particle charge, we are now in a posi- 
tion to sort them into categories which will enable us to choose those 
classes of events having the best resolution. A schematic version of the 
classification scheme is shown in Figure 3.5. The same scheme is used for 
both the high cutoff and the high ZC/ZI data separately. Only particles 
having a Cerenkov signal and at least one ion chamber signal were used 
in this analysis. All others were discarded. The selection limits were ori- 
ginally developed for use in other charge regions and minor "tweaking'' 
of the limits, although investigated, was deemed unnecessary because 
the small statistics of the analysis region precluded observation of all 
but large effects. The figures presented here serve as a diagnostic tool 
to ensure, for this first study, that the selections are reasonable, 
although not necessarily "optimal". 

The layout of the classification scheme can be understood in terms 
of particle fragmentation within the instrument, and the different 
categories by their varying amount of contamination by interacted 


Figure 3.5 


Schematic diagram of the selection "tree" used to categorize 
events on the basis of interactions. 



particles. If a particle penetrates the detector without any fragmenta- 
tion occurring, one would expect, at these high energies where slowing 
down is not important, that the charge estimates from the ion chamber 
or chambers, would match that of the Cerenkov counter. However, if the 
particle undergoes a nuclear interaction, we would expect the fragments 
to give a smaller signal than the original nucleus. This is a result of the 
Z 2 response of the detector. At the energies considered here the frag- 
ments usually have the same P and trajectory as the incident nucleus 
(e.g. Greiner et al. 1975), therefore only the charge dependence of the 
instrument response is involved in a comparison between the incident 
nucleus and its fragments. If 


Zmc = Zi + Z 2 (3.5) 

then 

Zinc 2 > Z a 2 + Z 2 2 (3.6) 

where is the intact incident nucleus charge and Z 1 and Z 2 are its 
fragments. Thus, by examining the separate estimates of particle 
charge, we can eliminate particles which have undergone a charge 
changing interaction whilt. passing through the instrument. Those 
events having the greatest number of consistent charge estimates 
should also include the least number of interactions and have the best 
resolution. 


The first major classification cf events in the "tree" of Figure 3.5 is 
on the basis of whether or not there are valid signals from ion chambers 
on both sides of the Cerenkov, i.e. two-module or one-module events. 
Having made this broad classification of particles we can now begin to 
eliminate interactions by examining the one module and the two module 
events separately. Looking first at the two module events, the first 
check lor interactions is accomplished by comparing the ion chamber 
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signals on both sides of the Cerenkov box. Between the two ion chamber 
modules there are the two aluminum honeycomb lids, or windows, each 
mainly consisting of ~1.24 g cm' 2 of aluminum, and the two Pilot 425 
Cerekov radiators with a mass of ~0.58fl g cm -2 of Lucite each (Binns 
I960). This is to be compared with the average interaction lengths for 
Sn-Ce of ~15 g cm' 2 in aluminum and ~6.0 g cm -2 in Lucite. Figure 3.6 is 
a histogram of 


21, - Zip 

TWDZ = 1 ZC 8 (3.7) 

for the high cutoff events passing the consistency checks explained in 
section 3.3, where we have used ZIj and ZI 2 to designate the two esti- 
mates of the charge derived from the mean of the ion chamber signals in 
each module of the sort described in 3.1 (7) (although for 90% of the 
events used here ZIj and ZI 2 are just the square root of the signal per 
unit pathlength) and ZC for the Cerenkov charge estimate. The histo- 
gram is centered at zero, showing that the ion chambers agree on the 
average. The ful width at half maximum is approximately 0.07 and is 
dominated by the ion chamber resolution with some contributions from 
the charge changing interactions present. If we assume that all the 
width is due to the ion chamber resolution, the TWDZ width implies that 
the ion chamber rms charge resolution is about 2%. (An ion chamber 
charge resolution of 1 charge unit at Z=40 roughly agrees with the 
observed distributions). The lower plot displays TWDZ as a function of ZC 
along with the approximate agreement criteria used. The quantity which 
was actually used in the selection was TWDZ + 0.5/ZC. The selection 
required it to be less than or equal to 0.06. 

It is still possible for interactions to have passed the TWDZ ion 
chamber agreement requirement. For example, in the charge 50-50 
region, a particle could have produced a fragment of 2 charge units or 
less and still be accepted by this cut. For these events, we would like to 
distinguish between two possible cases. First, the particle may have 


Figure 3.8 


Histogram shoving the effect of the two module agreement 
terion on the high charge, high rigidity data set. 

TWDZ = -- 1 ~ ZI 2 
ZC 

The erosaplot shows the selection as a function of ZC. 
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undergone the nuclear interaction after its passage through the Ceren- 
kov radiator. In this case ZC is still a valid charge estimator. The other 
possibility is that the interaction may have occurred before, or during, 
its passage through the Cerenkov, in which case ZC will not be a measure 
of the incident particle charge. In order to distinguish between these 
two possibilities a further selection is riade using 

ZDEL = — (3.6) 

where ZI^x = max [ ZIj, ZI 2 j. A histogram of this quantity for the high 
cutoff events which have agreement between modules is shown, in Figure 
3.7 along with a erossplot showing its dependence on ZC. The distribu- 
tion peaks at G.035 rather than zero both because of the normalization 
of ZI to minimum ionizing particles and because the use of Zl max biases 
the distribution towards higher values. The 8 GV requirement selects 
events which, because of the relativistic rise in the ion chamber 
response, have a larger ZI than that resulting from a minimum ionizing 
particle, thus causing the ratio of ZI/ZC to be greater than 1.0 for the 
events plotted. 

To understand the asymmetry of this distribution, we note that if a 
charge changing interaction of the type discussed above occurs before 
or within the Cerenkov the result will be a larger than ZC and thus 
a larger ZDEL than normal. On the other hand, an error in ZDEL on the 
low side is less likely to result from an interaction of the type considered 
here because it would imply that the Cerenkov signal is greater than 
those of the ion chambers on either side of it. In consequence, since we 
would like to eliminate interactions which will degrade ZC, the lower limit 
on acceptance is not crucial as long as the extreme outliers are elim- 
inated. Its value is -0.06. Using the uncertainty in the ion chamber ZI 
derived from the TWDZ distribution in Figure 3.6 and an estimated 



Figure 3.7 


Histogram of the agreement between Cerenkov and ion chamber 
charge estimates for high rigidity particles with two consistent 
ZI estimates. 


ZDEL = 


- ZC 
ZC 


The erossplot shows ZDEL as a function of ZC. 





Cerenkov resolution of 0.6% based on the widths of the individual tube ts 
average ratios (map corrected versions of Figure 3.1) we find that the 
expected ZD EL width should be ~0.02, consistent with what is shown. 
The upper limit of 0.06 is thus a selection at the ~1.25 sigma level 
(remembering that the peak is at 0.035). Using this value, we would 
expect some 10% of the "good", i.e. non-fragmented, events to be 
rejected (Category 6). Since Category 6 actually contains 20% of the 
TWDZ agreement events, half should be interactions. Adding this to the 
previous number outside the TWDZ limits (Category 3 and 4) implies a 
2fl% interaction rate overall, to be compared to the 30% expected from a 
simple interaction model which assumes a 17 g cm -2 interaction length 
in 5.6 g cm' 2 of aluminum traversed at a 20“ incident angle. The value of 
5.6 g cnf 2 includes not only the aluminum in the lid but also the "alumi- 
num equivalent" of the Cerenkov radiator which accounts for its shorter 
interaction length. Use of a 17 g em“ z interaction length instead of the 
15 g cm -2 more characteristic of the Sn-Ce region reflects the large 
number of particles in the distribution with charge near 40. Care should 
be taken in these comparisons since proton stripping reactions, which 
are still able to meet both the TWDZ and the ZDEL requirements, have 
not been excluded in the simple calculations done here. Those events 
accepted on the basis of both TWDZ and ZDEL (Category 5 in Figure 3 . 5 ) 
have three consistent charge estimates and should have the best charge 
resolution. The events in Category 6 on the other hand should be pri- 
marily those events for which ZC is measuring the combined fragment 
charges. 

We can also use the ZDEL parameter to further analyze those events 
which were rejected on the basis of TWDZ. Although these particles have 
undergone charge changing interactions between the two ion chamber 
modules, the events which interacted after passage through the Ceren- 
kov radiator may still have a valid ZC charge estimate. We would expect, 


since both the thickness of a module lid and a single radiator sheet are 
approximately the same fraction of an interaction length ( ~ 8%) that 
25% of these rejected TWDZ events ■would have a valid ZC, i.e. interacted 
in the lid after passage through both radiators and the first lid. The 
ZD EL histogram and cross-plot for the high cutoff events which failed the 
TWDZ test are shown in Figure 3.8. The histogram is similar to that 
shown in Figure 3.7 in its asymmetry. Note, however, that this distribu- 
tion peaks at a higher ZDEL value, reflecting the larger fraction cf 
events in this subset which interacted before or within the Cerenkov 
chamber. The same ZDEL cut used for Figure 3.7 was employed here. 
Again the lower limit is relatively unimportant. Use of 0.06 tor the upper 
limit can be justified on the grounds of attempting to select out the 
non-fragmented component which should have a distribution similar to 
the peak in Figure 3.7. The actual fraction of events in Category 3 
(interaction after Cerenkov) Is approximately 40%, higher than the 25% 
expected, which indicates that this category includes some particles 
which have interacted before or within the Cerenkov. There is an addi- 
tional selection made on those events fair which the direction of motion 
is "known". If the module entered first has a lower ZI than the second ZI, 
it is rejected without consideration of its ZDEL. This is a minor effect 
eliminating only one particle in the Sn-Ce region. Since these Category 3 
events have only two consistent charge estimates, the resolution is not 
expected to be as good as those in the analogous two module agreement 
Category 5. 

The one module events cannot be examined for interactions on the 
basis of module agreement. However, the ZDEL test can be used to 
select events with a consistent ZI and ZC. The ZDEL histogram for the 
one module, high cutoff events is shown in Figure 3.9. The most obvious 
feature is the large number of events with an unusually low Z'D'EL. 
Further investigation reveals that all of these events except 4 are "one 
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Figure 3.8 

Histogram of the agreement between Cerenkov and ion chamber 
charge estimates (ZDEL) for high rigidity particles with two 
inconsistent ZI estimates and a crossplot showing the ZC depen- 
dence. 





Figure 3.9 


Histogram and crossplot versus ZC of the agreement between 
Cerenkov and ion chamber charge estimates (ZDEL) for high 
rigidity particles having only one ZI estimate. 




radiator" events, i.e. events which., on the basis of their trajectory infor- 
mation, have penetrated only one of the two Incite sheets forming the 
Cerenkov radiator (RADI or RAD2 of Figure 2.2). These events would be 
expected to be of lower quality for several reasons. First, because of the 
instrument geometry, events in this category would also be expected to 
be wide angle and have only two hodoseope planes determining their tra- 
jectory. In addition, both planes are most likely on the same side of the 
Cerenkov box and, as a result, the computed position of the particle to 
the Cerenkov box has a larger uncertainty because of the "lever arm" 
from the hodoseope to the radiator. (Referring to Figure 2.2 again, to 
general, if the particle penetrated RAD2 and not RADI, it must have also 
penetrated the X3-Z3 HODO and the X4-Z4 HODO because of the coin- 
cidence requirement since it could not have gone through the X2-Z2 
HODO or the Xl-Zl HODO without penetrating RADI). If the particle 
actually did penetrate two radiators, but was assigned to the one radia- 
tor category iue to trajectory error, we would expect its ZC to be too 
high (because of an erroneous pathlength correction to the Cerenkov) 
with a resulting low ZDEL. Additionally, the map used to correct for 
Cerenkov areal response is based on two radiator events only and may 
only be approximately accurate for one radiator events. The ZDEL dis- 
tribution for the one module, one radiator events only is shown in Figure 
3.10. Note that there are two peaks, the one at ZDEL = —0.12, which may 
be attributable to the misassignment discussed above, and another at 
the location of the main peak. As can be seen from the crossplot, most 
of these low ZDEL particles are present only at lower ZC values. 

Figure 3.11 shows the ZDEL plots for the two radiator events of Fig- 
ure 3.9. This distribution is wider, rcns resolution of ~0.026, than that 
shown to Figure 3.7 for the two module events accepted on the basis of 
TWDZ. We note that here the definition of ZDEL differs slightly to its phy- 
sical meaning from the one used to the other selections. For two module 



Figure 3.10 


Histogram and crossplot versus ZC of t.he agreement between 
Cerenkov and ion chamber charge estimates (ZDEL) for the one 
radiator events in Figure 3.9. 




Figure ail 


Histogram crossplot of ZD EL for only the two radiator 
events of Figure 3.9 
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events, the use of ZI^., results in the preferential selection of the ion 
chamber signal from the module containing the non-interacted incident 
particle. As a result the ZDEL distribution does not include any values 
which were calculated with ZC measuring the incident particle and ZI the 
fragments (ignoring fluctuations). For one module events however, 
there is only one set of ion chambers and the selection of ZIj^ is 
superfluous. Thus it is equally likely, given an interaction, that ZC is 
measuring the incident particle and ZI the fragments as it is for the 
opposite case. The consequences of this do not affect the logic behind 
the selections. If the ion chamber measures the fragments and the 
Cerenkov the incident particle, ZDEL is low because ZC > ZI. Note how- 
ever that these events should have a valid ZC. If, on the other hand, the 
fragments are in the Cerenkov instead, ZDEL is high, as before. Thus the 
upper limit is again the more important one for eliminating interactions. 
For events on the low side of the ZDEL distribution, ZC should still be a 
valid charge estimate. 

In order to simulate this distribution, we plotted ZDEL for the e^eni * 
in Figure 3.7 (two consistent ion chambers) using ZIj or ZI 2 instead of 
Zlmax (the choice between them for each event being made at random). 
The distribution resembles the main peak shown in Figure 3.11, being 
wider, a sigma of ~ 03, than the ZDEL histogram made using We 

can therefore attribute this extra width to the use of the "fragmented” 
ZI. For the results presented here, because a comprehensive theory of 
one radiator events was outside the scope of this analysis, we did not 
separate them from the two radiator events and therefore they may 
have functioned as a possibly lower resolution component of the one 
module events. However, examination of the one radiator events 
separately reveal no significant difference between them and the two 
rcidiator events for the particles in the Sn-Ce region. 


- 71 - 


The same ZDEL criterion was used for the one module ZDEL plot. Fig- 
ure 3.9, as was used ior the others. As discussed previously, events with 
a low ZDEL should still have a valid ZC charge estimate. Although this is 
not true in the ease of a misassignrneat of two radiator events to the one 
radiator category, the crossplot shows that this is only important at the 
lower charges. Since the upper side of the ZDEL distribution is not 
affected by including events for which ZI is measuring the fragments, the 
same 0.06 limit is still appropriate. Using these values results in ?d (± 6) 
% of the events from the Sn-Ce region being placed in the interaction 
category (Category 2). To estimate this consider a simple model of frag- 
mentation which assumes a fragmentation pathlength of 15 g cm' 8 in 2.8 
g em“ 2 of aluminum traversed at a 45° incident angle. This calculation 
implies an 21% interaction rate. We should note here that we cannot dis- 
tinguish, for these one module events, between interactions occurring 
within the Cerenkov and those occurring between the Cerenkcv and the 
ion chamber for events moving from Cerenkov to ion chamber. As a 
result, a larger amount of contamination may be present Lr this 
category than is present in the two module events. 

The same selection ''tree" can also be applied to the events selected 
on the basis of ZC/ZI. The only category having sufficient resolution for 
ii. fusion in this analysis was Category 5 of Figure 3.5, the two module 
events which passed all agreement tests. Shown in Figure 3,12 is the 
TWDZ histogram and crossplot for the two module events. The width and 
number of events classified as interactions is similar to that of the two 
module-high cutoff events (Figure 3.6), The ZDEL plot for the particles 
which met the TWDZ requirement are shown in Figure 3.13 along with the 
ZDEL agreement criteria. The sharp cutoff in ZDEL on the high side is an 
artifact of the selection that ZC/ZI must be greater than 0.964. The 
cutoff is not perfectly sharp because the mean of ZI 3 and ZI 2 was used 
for the ZC/ZI selection while Zlm^ was used in ZDEL. 


& 
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Figure 3.12 

Two module agreement criterion for high ZC/ZI events as a 
function of ZC (lower) and as a histogram of TWDZ (upper). 

I 


r 
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Flgure 3.13 

Histogram and erossplot versus ZC o f ZDEL, the agreement 
between ion chamber and Cerenkov charge estimates, for two 
module high ZC/ZI events. 






3.6. Final Data Set 

The final decision of which of the non-interaction categories (1, 3, 
and 5 in Figure 3.5) to select for analysis was based on the resolution 
exhibited by histograms made using a modulo 2 superposition of the 
data with assigned Cerenkov charge between 49.0 and 59.0. The modulo 
2 superposition consists of binning the data according to the difference 
between its assigned charge and the nearest even integer. This has the 
effect of adding together the even (and odd) element peaks, thereby 
increasing the statistical accuracy with which the resolution can be 
determined. The method is useful here because the odd charge element 
abundances, for any proposed nucleosynthesis model in this region, tend 
to be lower than the adjacent even charge element abundances. As a 
result, the even elements should dominate the modulo 2 histograms. 

The categories actually used in the analysis are shown in Figure 3.14 
along with their associated modulo 2 histograms. Table 3.1 gives a list- 
ing of the number of high charge events meeting the consistency 
requirements of section 3.2 in each category, or bin of Figure 3.5. Also 
shown are the numbers obtained considering only the Sn-Ce region 
(49.0 ^ ZC :S 59.0). High cutoff and high ZC/ZI events are listed 
separately. In retrospect, the categories which were selected are not a 
surprise. The selected categories are those having the most information 
about each particle. The high cutoff, two module events have three con- 
sistent measurements of the charge and an independent estimate cf the 
particle energy. A Gaussian fit to their modulo 2 histogram, taking into 
account the spillover from neighboring elements, indicates a rms charge 
resolution of 0.47 ± 0.06 charge units, where the error is approximated 
using the formula for the uncertainty in the sigma of a gaussian distri- 
bution, i.e. -7 ^ — " with N the number of events in the distribution. The 
v2 N 

high cutoff, one module events, although having a modulo 2 peak which 
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Figure 3.14 


Final categories of Figure 3.5 selected for use in this data set. 
Histograms are plotted in 0.5 charge unit bins. Also shown are 
the modulo 2 histograms (in 0.25 charge unit bins) used in the 
selection. 

























is less regular and pronounced, still have a nominal resolution of 0.54 ± 
0.06. The two module high cutoff events in bin 3 which were salvaged 
from the TWDZ interaction events have the poorest statistics and cannot 
be shown to be statistically different from a uniform distribution. 
Nevertheless, there is a hint of resolution, a nominal 0.43 ± 0.12, and 
they might be expected to have at least some valid ZC charge estimates. 
The high ZC/ZI events, lacking the independent measure of energy, do 
not exhibit as good a resolution as their high cutoff counterparts having 
only 0.80 ± 0.05 charge unit resolution (cf. the two ''Bin 5" modulo 2 his- 
tograms in Figure 3.14). The lower quality categories not selected for 
use, 1 and 3 for the high ZC/ZI events, have a combined modulo 2 histo- 
gram containing 106 events which is consistent with a uniform distribu- 
tion containing the same number of events, the difference being 
significant at only the 26% level. 

The individual categories selected were normalized to iron by apply- 
ing the same selections to several days of iron data and then using a 
normalization factor to position the iron peak at 26. This same noimali- 
zation factor was then used to correct ZC for the high charge particles. 
This method worked we! for a! the high cutoff subsets giving normaliza- 
tions of 1.00142, 1.00361, and 1.00483 for Bins 1, 3, and 5 respectively. 
However, for the high ZC/ZI subsets, the normalization obtained in this 
way was incorrect. For the subset of these events actually used, those 
with two modules in agreement, the normalization was varied to obtain 
an approximately centered modulo 2 histogram. The final normalization 
factor was 1.0000, as compared to the iron result of 1.01521. The 
discrepancy between the iron and the modulo 2 normalizations for the 
high ZC/ZI subset has not been fully explained. However, it may be due 
to spectral differences, non~Z 2 effects in the ion chambers, or a combi- 
nation of both. Such effects, although not Lnsignific mt from iron (Z-26) 
to tin (Z=50), should be unimportant in the relative abundances of the 


elements from tin to cerium (Z=58). 

The final data set is shown in Figure 3.15. The main histogram, in 
G.25 charge unit bins, exhibits peaks at charge 38 and 40 serving to 
establish our charge scale. Possible charge dependent biases in the 
selections are such that relative abundances of widely separated 
charges should not be inferred from this plot. Inset a is an enlargement 
of the region of interest in 0.5 charge unit bins. One can clearly see 
peaks at the even elements soSn, 52 Te, 54 X 6 , sgBa, and sgCe. Also shown, 
as inset b, is the modulo 2 histogram, of the final data set. A Gaussian fit 
to this histogram, taking into account the spillover from neighboring 
elements, indicates a rms charge resolution of 0.55 ± 0.03, with the 
errors again estimated assuming a gaussian distribution. This resolution 
would be expected to be the result of similar contributions from pho- 
toelectron statistics (Garrard 1980) and residual mapping variations. 
The difference between it and a uniform distribution with the same 
number of events is significant at the greater than 99% level. The data 
set presented here has several differences from the one presented in 
Binns et aL. (1983). As the result of a reanalysis of the data, some addi- 
tional time periods were included (~4% increase) and a restriction on 
position in the Cerenkov radiator was relaxed (~12%). Additionally, all of 
the questionable trajectory events (58 in the Sn-Ce region) previously 
included in the data set, were eliminated here. These differences did not 
significantly change any of the relative abundances. 

The presence of a peak in the modulo 2 histogram, coupled with the 
presence of a peak at charge 3B, shows that our assumption of Z- 
squared scaling of the Cerenkov signal cannot be significantly in error. 
Since elements with even charge are more abundant than their neigh- 
boring odd elements, we are not likely to have an error of one charge 
unit, but an error of approximately two charge units cannot be immedi- 
ately ruled out. However, if there were a two charge unit error for any 
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Figure 3.15 

Histogram of the data from charge 35.0 to 60.0 in 0.25 charge 
unit bins. Although negligible over the limited range from 5 oSn 
to saCe, charge dependent biases in the consistency and energy 
selections may affect the relative abundances of widely 
separated charges. Insets show (a) the region of interest in 0.5 
charge unit bins and (6) a modulo 2 histogram of the data from 
charge 49.0 to 59.0 in 0.25 charge unit bins. 





one of the peaks in inset a, i.e. the Sn-Ce region,, while the peak at 38 
was correct, then the separation between even element peaks in this 
interval would be approximately 1.7 or 2.3 charge units (depending on 
the direction of the error) rather than 2.0, with the result that the 
modulo 2 histogram would not have such a well defined structure. Furth- 
ermore, such a large error in the charge estimate is not consistent with 
calculations of non-Z 2 Cerenkov effects by Derrickson et al. (1981). Ear- 
lier results from the Heavy Nuclei Experiment, using a different data 
subset have also shown the approximate validity of Z- squared scaling up 
to charge 40 (Binns et al. 1981b). 

Table 3.2 gives the abundances of ^Sh, 54X6, ^Ba, and 5B Ce normal- 
ized to our best estimate of the abundance of 5Z Te. These values have 
been obtained from fitting both even and odd elements, in the range 
from Z = 45.0 to Z = 60.0, to the 0.25 charge unit histogram using a 
Gaussian resolution function. The standard deviation of the assumed 
Gaussian was parametrically varied to obtain the best fit to the data. In 
addition, the abundances were constrained to be non-negative. The 
"quality of fit" was determined both by using minimum x 2 methods, 
which assume a Gaussian distribution for the uncertainty in each histo- 
gram bin, and by maximum likelihood methods, which assume a Poisson 
distribution. Results for the relative abundances of the even elements 
did not differ significantly between the two methods. The table values 
are derived from the v 2 fit using a value of 0.55 for the Gaussian stan- 
dard deviation. The uncertainties are the limits, for the given parame- 
ter, at which the x 2 can be made equal to the minimum x 2 plus one by 
allowing the other parameters to vary (see e.g. Bevington 1969). This 
should correspond, approximately, to the one standard deviation errors. 
The "best fit" abundances have a x 2 of 40.47 for 43 degrees of freedom. 
There is a 72.2% probability of obtaining a x 2 this high or higher by 
chance. The best fit value of the charge resolution has x 2 + 1 limits 



TABLE 3.3 


Element 

Abundance 
(Fit Results) 

Abundance 

( gaTe = 1) 

so Sn 

41.23 ± 9.47 

1.65 ± 0.30 

52Te 

24.93 ± 6.90 

1.00 ± 0.26 


25.53 ± 6.73 

1 .02 ± 0.27 


44.91 ± 7.40 

1.80 ± 0.30 

ss^e 

20.27 ± 6.23 

0.01 ± 0.25 











extending from approximately 0.44 to 0.80 change units. Also indicated 
in Table 3.2 are the actual abundance values derived from the x 2 fit. 

No corrections have been applied to the data for fragmentation in 
the instrument. Employing a simple model of fragmentation, which uses 
an empirical geometrical cross section, a, for fragmentation of the 
incident nuclei of the form 


a — n rf [a^ / 3 + A^ /3 — bj (3.9) 

where r Q = 1.35 fm, b — 0.83, (Westfall et al., 1979) and At and Ag are the 
mass numbers of the target (aluminum) and incident nucleus respec- 
tively, and which assumes an equal probability for production of all 
lower charge fragments, the relative adjustment factors for the the even 
charge nuclei in this limited charge range are less than 5%. 



TV. DISCUSSION 


4.1, Overview 

The measurements presented in the previous chapter are the first in 
this charge range to exhibit resolution of the even charge elements. In 
this chapter we will discuss these results in light of the nucleosynthesis 
processes believed responsible for the formation of the ultra-heavy 
cosmic rays. The first three sections will enumerate the complications 
and uncertainties which arise In comparing the measured cosmic ray 
flux to the theoretical results of nucleosynthesis. The next two sections 
will compare our results to the various models and to previous measure- 
ments in the same charge interval. Finally, we will conclude with a brier 
discussion and summary of our findings. 

4.2. Nucleosynthesis Source Models 

The first obstacle to a comparison between cosmic ray measure- 
ments and nucleosynthesis theory is determining what the theoretical 
results should be. Specifically, what are the abundances which are the 
result of "normal'' nucleosynthesis? The canonical reference to which 
cosmic ray abundances have been compared for years has been one of 
the set of periodically updated versions of the "solar system" abun- 
dances of Cameron (I960, 1973, 1902b). These abundances are based 
mainly on measurements made on type Cl carbonaceous chondrites, 
with solar abundance data and some nucleosynthesis theory being used 
to fill in the gaps. Carbonaceous chondrite meteorites are chosen 
because they are believed to be the most representative of the primitive 
solar nebula. Recently, a new abundance compilation was published by 
Anders and Ebihara (1902) which is also based on type Cl carbonaceous 
chondrites and has a larger data base (a factor of two or more) than 


that of the most recent Cameron tables. In addition, the new compila- 
tion includes estimates of the errors on the assigned abundances, a 
feature which has been lacking previously. 

Figure 4.1 is a logarithmic plot of the even charge element abun- 
dances in the Sn-Ce region using the compilation of Anders and Ebihara 
along with their estimated errors. Also included on the plot are the 
values of Cameron (1982b). As can be seen, the agreement between the 
two compilations is good except for an approximately 24% decrease in 
the Anders and Ebihara values for tellurium (agTe) and xenon ( 54 X 6 ) 
compared with those of Cameron. The value for Te has been decreased 
because of a systematic error in the data upon which the old value was 
based (see Anders and Ebihara 1982 and references therein). The 
decrease in Xe follows from the Te decrease because the abundance of 
Xe, a noble gas, is not based on actual meteorite measurements, which 
exhibit a great deal of variability from meteorite to meteorite, but 
rather on the results of a fit to the "Te-I-Cs-Ba peak". Cameron also uses 
a similar type of interpolation to determine hisXe abundance. 

Besides the solar system abundances, we are also interested in a 
comparison with the results of r- and s-proeess nucleosynthesis. In 
order to separate the contributions of the s- and the r-proeess in the 
solar abundances, we must decompose them according to the scheme 
discussed in chapter I. Although there have been several recent decom- 
positions (e.g. Israel et ai. 1981, Blake and Margolis 1981, Cameron 
1982a), we have chosen to use the results of Kappeler et ai. (1982) as our 
model s-process. This calculation uses new and improved measurements 
of neutron capture cross sections of important s-only and magic number 
nuclei along the s-process path. In particular. l 3 e Ba, with a magic 
number of neutrons, is included in this group. 

The crN curve of K&ppeler et ai. is the result of an s-process 
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Figure 4.1 

A his tog ram of the abundances of the even charge elements in 
the goSn to aaCe region. The vertical axis is logarithmic with the 
abundance of Si defined to be 10 8 . The solid histogram uses the 
results of Anders and Ebihara (1982) while the dashed lines are 
the results of Cameron ( 1982b). The error bars are those of 
Anders and Ebihara. 
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calculation which assumes a steady neutron flux and. a two component 
exponential distribution of exposures. The input parameters, consisting 
of the fraction of iron seed nuclei exposed to the two components and 
the two mean values of neutron exposure, were varied to obtain a best 
fit to the Cameron 1982b abundances of selected s-only nuclei. Although 
i the model was fit to the the Cameron abundances, we can still use this 

s-proeess in conjunction with the new Anders and Ebihara abundances. 
Figure 4.2 is a plot of the Kfippeler et al. crN curve as a function of mass 
number, adapted from their paper. Indicated on the plot are the empiri- 
cal values of crN for s-only and predominantly s nuclei which illustrates 
how well the curve fits the data. The open circles are the abundances of 
Cameron (1982b) multiplied by the neutron capture cross sections 
reported in the K&ppeler et al. paper. The filled circles use the same 
neutron capture cross sections with the eb’indanf'es of Anders and 
Ebihara. The error bars attached to the Cameron points use only the 
uncertainties in the cross sections (Kdppeler et al. 1982) and do not 
include any abundance uncertainties. As a result similar error bars 
would also apply for the Anders and Ebihara points:. Qualitatively, the 
curve appears to fit both sets of data points equally well with no gross 
differences. In particular, the fit in the Sn-Ba region is improved using 
the Anders and Ebihara abundances as a resuit of the decrease in th-e 
value for Te, affecting the s-only isotopes 1Z2 Te, 123 Te, and 124 Te. (This 
decrease was actually recommended by K&ppeler et al. in their original 
paper.) We note in passing that the "ledge-precipice” structure is a 
characteristic resulting from the "bottlenecks" at the magic number 
nuclei with their extremely small neutron capture cross sections. 

Figure 4.3 shows the s- and r-process values in the Sn-Ce region 
which result from subtracting the Kappeler et al. s-proeess from the 
Anders and Ebihara abundances on an isotope by isotope basis. Also 
shown for comparison is the r-proeess which results if the original 
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Figure 4.2 


A plat of the crN curve of Kappeler et a l. (1982) along with, the 
empirical oN product for s-only and predominantly s isotopes. 
The values of Cameron (1982b) are shown as the open circles 
with error bars which include only the neutron capture cross 
section uncertainties. The filled circles are the abundances of 
Anders ard Ebihara (1982) and would have the same uncertain- 
ties. The neutron capture cross section, a, is measured in miM- 
barns and N, the s-process abundance, is based on the usual Si 
= 10 e scale. 






figure 4.3 


The s-proeess (left) and r-process (right) components of the 
solar system in the 30 Sn - 3B Ce region using the Kappeler et a L. 
s-process and subtracting it from Cameron (1982b) (dashed 
lines) and the Anders and Ebihara (1982) (solid lines) solar sys- 
tem values. The errors bars are estimates based on the uncer- 
tainties quoted in Anders and Ebihara. The p-process nuclei 
have been excluded in this analysis. 








Cameron (1982b) solar system values are used (from K&ppeler et al.). 
Table 4.1 shows the details of the decomposition used here for the even 
charge elements between soSn and eoNd. Neodymium is included because 
it was used in calculating uncertainties later on (see section 4.4). The 
isotopes which could only be made in the p-process were not included in 
either the r- or the s-process. As a result, the sum of the r-process and 
the s-proeess does not necessarily equal the solar system value. Addi- 
tionally, for the s-only isotopes, the value derived from the crN curve 
(Table 7 of K&ppeler et al.), and not the actual isotopic abundance, was 
used. R-proeess abundances of these s-only nuclei were defined to be 
zero. The largest difference between the solar system and the sum of r- 
and s-process elemental abundances resulting from this procedure is 5% 
for Sn with the other elements having less than a 1% effect present. 

Table 4.2 summarizes the elemental abundances. The uncertainties 
in the solar system abundances are those of Anders and Ebihara. For 
the s-process, in order to test the sensitivity of our results to the form 
of the s-process used, we have estimated the limits of variation probable 
by using the percentage uncertainty in the neutron capture cross sec- 
tion for a given isotope as the percentage uncertainty in the s-process 
abundance of that isotope. This is reason ■’hie on the grounds that the s 
abundance of the nuclide with mass number A is determined by dividing 
the oN curve, which is relatively invariant, by the cross section, u k . A 
more appropriate value to use, in light of the equation for the cN curve 
for an exponential exposure of neutrons is the percentage uncertainty 
in 




<j a t 


(4-1) 


where r is the mean neutron exposure and is assumed to have no uncer- 
tainty here. For an exponential distribution of exposures, the 


Table 4.1 



Table 4.1 - 50 Sn 


A 

Solar 

System 

... 

s- 

proeess 

r- 

process 

% 

error 

112 

0.0386 

- 


- 

114 

0.0256 

- 

- 

- 

115 

0.0145 


- 

- 

116 

0.565 

0.458 

- 

20 a 

117 

0.296 

0.105 

-0. 191 

7 b 

110 

0.929 

0.655 

0.274 

8 b 

119 

0.329 

0.156 

0.173 

15 b 

120 

1.24 

0.749 

0.49 

30 b 

122 

0.174 

- 

0.174 

- 

124 

0.215 


0.215 


Total 

3.82 

2.123 

1.52 

- 


Percentage uncertainties apply to the s-process abundances. 
11 Kappeler et aL. (1902) 
b Alien, Gibbons, and Maekiin (1971) 
e Uncertainty of 30% adopted (50% for Xe) - see text. 
d Uncertainty in (4.1) used with 
r = 0.092 mb -1 
Oia8 3a = 4.22±0.25mb 
(7i« Ce = 11.5±0.6mb 
(Kappeler et a L. 1982). 















Table 4.1 - ^Xe 



124 

126 

128 

129 

130 

131 

132 
134 
136 

Total 


Solar 

System 

s - 

process 

r- 

process 

0.00496 

- 

- 

0.00483 

- 

- 

0.0939 

0.108 

- 

1.20 

0.0603 

1.14 

0.189 

0.176 

- 

0.941 

0.0619 

0.879 

i.15 

0.254 

0.90 

0.421 

- 

0.421 

0.34 

- 

0.34 

4.35 

0.660 

3.68 


% 

error 


50 a 
50 c 
30 a 
50 c 
50 c 




































Table 4.1 - M Nd 


A 

Solar 

System 



s- 

proeess 

r- 

process 

% 

error 

142 

0.227 

0.151 

- 

19 a 

143 

0.101 

0.0292 

0.072 

30 0 

144 

0.199 

0.112 

0.087 

30 c 

145 

0.0694 

0.0158 

0.0536 

30 c 

146 

0.144 

0.0607 

0.083 

30 c 

148 

0.0477 

- 

0.0477 

- 

150 

n na«R 

- 

Q.0468 

- 

total 

0.835 

0.369 

0.390 

- 











T 


Table 4.2 

Source A bundances (Si s 10 6 ) 


Element Solar System s-process r-process 


3.82 

± 

9.4% 

4.91 

± 

12% 

4.35 

± 

13% 

4.36 

± 

4.5% 

1.16 

± 

5.178 

0.835 

± 

7.0% 


2.123 ± 12% 

| 1.52 ± 25% 

0.920 x 6.4% 

3.96 ± 14% 

0.660 ± 23% 

3.68 ± 15% 

4.11 ± 3.8 % 

0.29 ± 78% 

0.765 ± 2.7% 

0.390 ± 16% 

0.369 ± 13% 

0.390 ± 15% 


o 








abundance of the nuclide with mass number A along the s-process path 
is proportional to 


1 ^ 

— n 

a k i«S6 


1 + 


a, r 


-l 


(4-2) 


with tTj the neutron capture cross section of the nuclide with mass 
number i along the s-process path. (e.g. see K&ppeler et cd. 1982). Since 
we are interested only in the uncertainties in the nuclide's relative 
abundance, we ignore the uncertainty from all terms in the product 
which are common among the nuclei under consideration leaving only 
the term (4.1). The use of (4.1) is only important for the magic number 
nuclei where o a t <1. We have calculated the uncertainties in the iso- 
topes I38 Ba and I40 Ce using (4.1). For all the others we used the uncer- 
tainty in a*, alone. The isotopic uncertainties assumed are listed in Table 
4.1 in the last column. The values are of three types. For those listed in 
Kappeler et a!., we used their uncertainty. For those uncertainties not 
listed in Kappeler et at. we used the uncertainties listed in Allen et a I . 
(1971) since, for a number of the nuclei, the cross sections were equal to 
the values inferred from the <rN - N table in Kappeler et cd. For those 
cross sections where neither compilation listed uncertainties we used 
the minimum value of 30% suggested by Kappeler et cd. except lor Xe 
where 50% seemed more appropriate on the basis of the listed uncer- 
tainties. The uncertainties in the r-process are just the uncertainties in 
the s-process and solar system values added in quadrature for the iso- 
topes involved. 
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4.3. Preferential Acceleration 

A further complication to the determination of the cosmic ray 
source abundances arises from possible preferential acceleration 
effects. Comparisons of the derived cosmic ray source abundances with 
solar system abundances show indications of a possible first ionization 
potential effect of the type discussed in Cassd and Goret (19?0) and 
references therein. Those elements having a low first ionization poten- 
tial also have enhancements in the ratio of cosmic ray source abun- 
dance to solar system abundance which may be indicative of the pre- 
ferential acceleration of those elements easiest to ionize. 


Figure 4.4 is a plot of the ratio of the cosmic ray source abundance 
to the solar abundance as a function of the first ionization potential as 
taken from Brewster, Freier, and Waddington (1983). The ratio is defined 
to be 1 for Fe. It can be seen that there is a definite, but not perfect, 
correlation between the twc quantities. The work of Cassd and Goret 
was based only on results in the charge region below iron. Nevertheless, 
results in the charge 36 to 40 region (Binns et ai. 1903) have shown that 
this correlation seems to be present for the higher charges also. Both 
sets of data are included in the figure. 

The actual form of the first ionization potential dependence has not 
been fully established. Some models involve a step function with the 
effect being discontinuous at a value of ~ 9 eV, Others use an exponen- 
tial dependence. But even here the numerical values used in the func- 
tion depend on what range of ionization potentials is used for the fit. In 
particular, the two exponentials indicated on the figure by the straight 
lines are results obtained both with and without including the high first 
ionization potential elements He and Ne. For our comparison, we have 
used the "FEPl'' form of Brewster et al. 


R = 9.3 1 exp' 


(4.3) 


Figure 4.4 


Plot of the ratio of the cosmic ray source abundance to the 
solar system abundance as a function of the first ionization 
potential of the element. Also shown are two possible fits to the 
trend in the data. From Brewster, Freier, and Waddington 
(1903). 



IONISATION POTENTIAL 




which was chosen to fit the data with first ionization potential, i, 
between 5 and 15 eV. R is the ratio of cosmic ray source to solar system. 
Table 4.3 shows, for the even elements in the Sn-Ce region, the values of 
R given by this fit along with the first ionization potentials assumed. In 
order to assess the sensitivity of our results to this effect, we have 
adopted an uncertainty of 15% for FIP adjustments. 

4.4. Cosmic Ray Propagation 

Another important effect which alters the abundances observed at 
earth is due to cosmic ray propagation. As the cosmic rays pass through 
the interstellar medium, they undergo nuclear interactions with the H 
and Ke of which it is primarily composed. These interactions have the 
effect of reducing the flux of the more abundant nuclei by spallation and 
of increasing the flux of the rarer nuclei with the fragments of the more 
abundant ones. For this study, we use the propagation calculations of 
Brewster, Freier, and Waddington (1983). These calculations model the 
cosmic ray propagation effects using matrix methods and the "leaky 
box" formalism of CowsBc et a l. (1967). In the steady state leaky box 
model, the observed flux of cosmic ray species i is related to its rate of 
production at the source by the continuity equation (e.g. Stone and 
Wiedenbeck 1979) 


0 = 4 + St 1 - h (4.4) 

j "ij 

where $j is the observed flux of species i, q, is its rate of production at 
the source, Aj is the total interaction length of species i, A :J is the 
interaction length for production of species I from species j, and the sum 
over j is over all species able to fragment into i. We nave ignored energy 
losses in (4.4). The model is essentially an equilibrium calculation in 
which nuclear fragmentation and decay, along with escape from the 


Table 4.3 

FIP Values 

Element 

I(eV) 

R(FTF1) 

goSn 

7.34 

1.12 

52^e 

9.01 

0.695 

64Xe 

12.13 

0.283 

M Ba 

5.21 

P.08 

saCe 

5.60 

1.06 

floNd 

5.50 

1.91 


Ionization potentials and FTP 1 from Brewster (private communi- 
cation and Brewster, Freier, and Waddington 1983). 
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galaxy, deplete the source nuclei, which are simultaneously being 
replenished both by fragmentation of nig her mass nuclei and by a con- 
tinuously operating source. Since it is a steady state model, qj and 
are respectively proportional to the abundances of species i at the 
source and as observed at earth. For details of the Brewster, Freier, and 
Waddington calculation one is referred to their article. A few points 
deserve mention here, however. First, no account is taken of energy loss 
during transport in their calculation, which should be valid for the ener- 
gies of 2 GeV/nueleon and higher that we are working with. Second, the 
cosmic rays are considered to have a mean escape length of 5.5 g cm" 2 
in an interstellar medium composed of pure H. This results in an 
exponential distribution of pathlengths with this same mean. Finally, 
although a complete propagation involves isotopes and not elements, 
because of the large number of possible nuclides to consider, this calcu- 
lation uses interaction cross sections weighted according to the 
assumed source isotopic composition and propagates elements only. 

Table 4.4 lasts the results of propagating the various sources dis- 
cussed in section 4.1 (Brewster, private communication). The propaga- 
tion included all the elements with Z& 50 (see Appendix B). The values 
for the s-process are not the results of an actual propagation but were 
derived by subtracting the propagated reprocess elemental abundances 
from those of the propagated solar system Although this procedure is 
not exact, any errors made should be small compared to the uncertain- 
ties involved. The quoted error limits ore estimates based on the sou ee 
and propagation uncertainties (Table 4.2 source abundance uncertain- 
ties were used. FIP values include a 15% uncorrelated uncertainty in the 
FIP source adjustments). For calculating errors, a simple model of the 
form 



(4.5) 


M 


Table 4.4 

Propagated Fluxes (Arbitrary Units) 


NO FTP 

Element 

Solar System 

s- process 

r-process 

5oSn 

«Te 

5iXe 

56® a 

5 fiCe 

7.45 ± 16% (0.58) 

6.86 ± 13% (0.79) 

5.86 ± 13% (0.81) 

5.25 ± 7.7% (0.86) 

1.47 ± 8.9% (0.81) 

3.83 ± 15% (0.62) 

1.76 ±18% (0.57) 

1.54 ± 20% (0.45) 

4.505 ± 6.4% (0.94) 

0.882 ± 6.7% (0.89) 

3.62 ± 22% (0.47) 

4.90 ± 14% (0.B7) 

4.12 ± 16% (0.94) 

0.745 ± 38% (0.40) 

0.588 ± 16% (0.68) 


Element 

Solar System 

s-proeess 

reprocess 

5oSn 

8.02 ± 19% 

(0.60) 

4.64 ± 19% 

(0.57) 

3.38 ± 23% 

(0.57) 

saTe 

5.27 ± 18% 

(0.70) 

1.81 ± 26% 

(0.30) 

3.46 ± 19% 

(0.06) 


3.36 ± 23% 

(0.39) 

1.74 ± 31% 

(0.11) 

1.62 ± 31% 

(0.68) 

56® a 

10.4 ± 15% 

(0.90) 

9.16 ± 16% 

(0.96) 

1.24 ± 44% 

(0.50) 

5flCe 

2.56 ± 15% 

(0.86) 

1.584 ± 15% 

(0.92) 

0.976 ± 20% 

(0.76) 


Propagated values from Brewster, private communication. See also 
Brewster, Freier, and Wadding-ton (1983). The number in parentheses 
is the fraction of the observed abundance assumed to be surviving pri- 


oinn 
























was assumed. In (4.4), if has no contribution from the secondary com- 
ponent, the summation term, then we have 

♦j = Ajqj (4 6) 

and $j is all surviving primary, i.e. no secondaries. This is approximately 
true for Bi which has no abundant elements above it on the charge scale. 
For an element i for which j is the major secondary contributor, we can 
ignore all terms other than j in the secondary sum and substitute Ajqj for 
giving us (4.5). In making this approximation we thus ignore the con- 
tributions made to element i from ah other elements except the source, 
or primary, component of j. It was assumed that the total interaction 
cross sections could have an average error of 5% (Letaw e£ al. 1983) and 
that the partial cross sections had uncertainties of 30% (although they 
could be in error by as much as 50%, see e.g. Brewster, Freier, and Wad- 
dington, 1983). Additionally, in the calculation of the uncertainty con- 
tributed by the fragments from higher charges to the observed flux, 
approximated by the second term above, the uncertainty of the next 
higher charge even element was used under the assumption that it was 
the main contributor. The number in parentheses for each entry is the 
fraction of the observed flux attributed to the surviving primary 

4.5. Comparison of the Data with the Models 

Figure 4.5 is a comparison between the results of this study and the 
propagated even element fluxes in the Sn-Ce charge region. No first ion- 
ization potential adjustments have been applied to the sources and all 
the model distributions have been normalized so that the abundance of 
52 Te is defined to be one. Because we have constrained our analysis to a 
limited charge region, we are able to vary the normalization between th 0 
model and the data to obtain a 'best fit", defined here as a minimum in 
the x 2 , fo r the five elements under consideration, 5 oSn, 52 Te, 54 Xe, 5 eBu, 
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Figune 4.5 

Comparison of the data with the results of propagating a solar 
system type source, a pure s-proeess source, and a pure re- 
process source derived from the abundances of Anders and 
Ebihara (1982) and the s-proeess of Kappeler et al. (1982). 



A. \ 
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and gflCe. The adjustment of the normalization means that the x 2 has 
four degrees of freedom. The data are shown with this best fit normali- 
zation applied. Note that the vertical scale is logarithmic. The x 2 calcu- 
lation and the error bars indicated on the data do not include the 
source uncertainties, only those due to the data points, although, as will 
be seen later, this does not change the nature of the results presented 
here. 

The three panels show, from left to right, the results of propagating 
solar system, pure s-proeess, and pure r-process sources. Note that best 
agreement seems to be with the pure s-process source if no first ioniza- 
tion potential biases are included. Table 4.5 shows the x 2 values and the 
probability of obtaining a x 2 that high or higher by chance. 

Figure 4.6 is the same as figure 4.5 except that the sources have had 
first ionization potential (FEP) adjustments applied before propagation. 
In this case, the best fit is obtained for a solar system mixture of r- and 
s-process material with a x 2 of 3.99 for four degrees of freedom which 
corresponds to a 40.6% probability of occurring by chance. The values 
for the other sources are shown in Table 4.5 under FTP. 


Although the previous figures make a pure r-process source highly 
unlikely on the basis of our measurement, it does not rule out the possi- 
bility that a significant fraction of the cosmic ray source may be the 
result of r-process synthesis, in order to examine this question, we have 
looked at the x 2 as a function of the relative amounts of r- and s-process 
material present at the source by constructing the composite abun- 
dance distribution Xj. 


Xj = fSj + ( 1 ~ f)Rj (4.7) 

where Sj denotes the s-process abundance of element i, Rj denotes its r- 
process abundance, and f is a parameter which varies fromO to 1. If the 


r-proeess 


Table 4.5 
X&in Values 


Solar System 


NO FIP 39.59 (<0.00001%) 11,04 

+ERR0RS 38.45 (<0.00001%) 11.61 


FIP 34.37 
+ ERRORS 10.65 


s-process 

5.32 

(25.6%) 

5.24 

(26.4%) 

16.50 

( 0.36%) 

7.63 

(10.6%) 


"+ ERRORS" includes uncertainties in the source abundances 
and propagation uncertainties. Numbers in parentheses are the 
probability of obtaining a x 2 that high or higher by chance. 
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Figure 4.6 

Comparison of the data with the results of propagating a solar 
system type source, a pure s-process source, and a pure r- 
proeess source derived from the abundances of Anders and 
Ebihara (19B2) and the s-process of Kappeler et al. (1982) and 
adjusted for first ionization potential effects as discussed in the 
text. 
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Sj have the correct ratios to each other for the s-process and the Rj have 
the correct ratios for the r-process, then the resultant Xj represents a 
mixture whose relative elemental abundances vary from pure r-process 
to pure s-process as f goes from 0 to 1. hi particular, if we chose the S, 
and Rj of Table 4.4, then if f = 0.5 


3Q = 0.5 (Rj + Si) 


(4.6) 


and we have the solar system distribution to within a normalization fac- 
tor. Since we will only be concerned with the relative abundances of the 
elements in this region, the added factor of 0.5 is inconsequential. A 
parameter that is more physically meaningful than f is 


1 - f _ (r/s)^ 

i (r/sJss, 


(4.9) 


where the second equality, the ratio of r-process to s-process matenal in 
Xj compared with the same ratio in the solar system, follows if and S, 
are chosen from Table 4.4, i.e. so that SSj •- Rj + S* where SSj are the 
solar system abundances. When f = 0.5, 77 = 1 and the distribution con- 
tains the same fraction of s- and r- process material as does the solar 
system. If f = 0.4, the Xj distribution contains 77 = 1.5 times the *olar 
ratio of r- to s-process material. 

Although the above formalism will work for the source distributions, 
it is not immediately obvious than one can do the same for the pro- 
pagated fluxes. However, since the propagation operation is, in essence, 
a matrix multiplication, a linear combination of sources is the same 
linear combination of propagated fluxes. This is only true exactly if one 
is working with isotopes because the isotopic cross sections for a given 
element are not all equal and, since the isotopic as well as the elemental 
abundances vary, this means that the elemental cross section are a 
function of f in the Brewster, Freier, and Wad ding ten propagation and 


therefore so are the matrix elements. Nevertheless, the dependence on f 
is weak in the weighted cross sections used, generally showing variations 
of less than 5% from pure r-proeess to pure s-process material (Brew- 
ster, private communication). As a result, the errors should not be 
sig nificant here. 

Using the above formalism, the x 2 was found for each value of f 
between 0 and 1 in the same manner as described previously, i.e. using 
the best fit normalization. Figure 4.7 is the result for the ease where no 
"TP adjustments were applied. The pure r-process is at the left, the solar 
system at f = 0.5, and the pure s-process at the right. The left hand 
scale indicates the x 2 values with selected x 2 significance levels indicated 
by the labeled dashed lines. As can be seen from the figure, the best fit 
is for an f of 0.83 which corresponds to an r-process to s-process ratio of 
0.20 times that of the solar system. The x 2 + 1 uncertainties on this 
quantity are 0.725 to 0.945 in f, corresponding to an r- to s-process ratio 
which is from 5.5% to 38% of the solar system value. If we include the 
errors on the propagated sources listed in Table 4.4, the results are 
essentially unchanged (Table 4.5 "NO FIP + ERRORS”). The x 2 curve, to 
within the accuracy depicted, is the same for the range shown. The pure 
r-process x 2 drops from 39.6 to 38.4, s 'll a significant difference at the 
greater than 99% level. This relative insensitivity to the model uncer- 
tainties implies that the data errors dominate the calculation. The abso- 
lute uncertainty on the Sn point, for example, changes only from 0.38 to 
approximately 0.42 if the r-process, best fit errors are added in quadra- 
ture. The Ba uncertainty remains unchanged. 

We can do the same calculation for the FIP adjusted sources. These 
results are shown in Figure 4.8. The format is the same as in Figure 4.7 
with the solid curve indicating the x 2 value when only the data point 
errors are considered. The dashed curve results if the uncertainties in 
the propagation model are included (Table 4.4). This greater senadtiviiy 
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figure 4.7 

X Z of the fit between the data and a given mixture of r and s- 
proeess material as a function of the mixrng parameter, f. Solar 
system material corresponds to f = 0.5 while a pure r-proeess 
corresponds to f = 0 and a pure s-proeess to f = 1. 
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Figure 4.8 

X 2 of the fit between the data and a given mixture of r and s- 
proeess FIP adjusted material as a function of the mixing 
parameter, f. Solar system material corresponds to f = 0.5 while 
a pure r-process corresponds to f = 0 and a pure s-proeess to f 
= 1. The solid curve includes only the data uncertainties. The 
dashed curve includes uncertainties in the model as presented 
in Table 4.4. 



of the calculated x 2 to the uncertainties on the model is mainly the 
result of the FEP adjustment increasing the Ba abundance by a factor of 
approximately 2, resulting in larger absolute uncertainties in this quan- 
tity, Nevertheless, the general character of the curve does not change. 
The minimum in bath cases occurs at 0.40 which is an r-process enrich- 
ment of 1.50 over the solar system value. Again, the x 2 + 1 uncertainty 
limits are 0.28 to 0.55 in f or 2.6 to 0.82 in tj. (If we consider only the 
data errors, these limits become 0.30 to 0.51 for f and 2.3 to 0.96 in 77 ). 

The sensitivity of the calculation to the different element abun- 
dances is indicated by Figures 4.9 and 4.10. In these, the quantity 


Oj = 


Dj-X; 


(4.10) 


is plotted as a function of f. Here, X< is the r-s combination for element i 
with the best-fit normalization applied, Dj is the data value, and Cj is the 
uncertainty in Dj. For the illustrated curves, <7j includes only the data 
errors. Each curve is labeled with the appropriate element. Figure 4.9 
addresses the non-FIP case, while Figure 4.10 applies for the FIF 
adjusted values. Note that in both cases Ba shows the greatest degree of 
variation and that, additionally, a value of f can be found for which all 
elements are in reasonable agreement with the theory, as could be 
expected from the x 2 analysis. Cerium is the only element which 
remains significantly overabundant both with and without FIP. However, 
the effect is less than two sigma for all acceptable values of x 2 and is 
therefore marginal at best. 

Another, slightly different, way of viewing the data can be seen in 
Figure 4.1 1. In this diagram, the ratio of Ba to Te is plotted against the 
ratio of Sn to Te. Since Sn and Ba are both primarily s-process ele- 
ments, while Te is primarily r-process, this plot is effectively a display of 
the correlation between these two measurements of the amount of r- 
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Figure 4.9 

Contributions to the x 2 for each element as a function of the 
mixing parameter f discussed in the text. No first ionization 
potential adjustments were applied. 
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Figure 4.10 

Contributions to the x 2 for each element as a function of the 
mixing parameter f discussed in the text. Source abundances 
were altered according to first ionization potential biases. 



Figure 4.11 


The results of our measurement and of the propagation calcula- 
tion of Brewster, Freier and Wad ding ton (1983) applied to r, s, 
and solar system sources both with and without first ionization 
potential (FIP) effects included are shown. The solid and dashed 
contours illustrate the 68% and 50% significance levels based on 
a x 2 calculation using the three elements soSn, 52 Te, and sgBa. 
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and s-process material present. The data point is shown as the large 
filed circle surrounded by the 50% (dashed) and 60% (solid) x 2 error 
contours. Their elliptical shape is the result of correlations arising from 
using Te as the normalizer for both axes. Also shown on the plot are the 
results from the Brewster, Freier, and Waddington propagations. The 
dashed line connects those points resulting from the pure r-process, the 
solar system, and the pure s-process, with no first ionization potential 
adjustments applied. The solid does the same for the FIP sources. These 
lines are "mixing lines" with position along the line indicating the rela- 
tive amounts of r- and s-process material present, similar to the Xj used 
previously. This figure is consistent with the results obtained using all 5 
elements. We are in agreement with both a solar system source with FIP 
applied and with an s-dominated mixture if no FIP is used. Using only 
the three elements ^Sn, ogTe, and sgBa, we are inconsistent with a pure 
reprocess source at the greater than 99% level whether or not first ioni- 
zation potential effects are included. Consideration of model uncertain- 
ties reduces the significance of the difference with r-process plus FTP to 
the 90% level but leaves the no FIP r-process significance essentially 
unchanged. 

4.6. Comparison with Other Measurements 

Table 4.6 lists our results along with recent measurements from two 
other experiments, the Ariel VI electronic detector and a balloon borne 
plastic track detector (Fowler et al., 1901). The balloon values are ''sub- 
ject to significant and rapidly charge-dependent corrections due to 
threshold effects for the tracks in plastic" (Fowler et al., 1901). For the 
Ariel VI data, no claim of resolving the individual even elements in this 
charge range has been made. 

Although the balloon data are consistent with a uniform distribution 
in this limited charge region, the Ariel VI results do posses a significant 


Element 


HEAO 


3le 4.8 


mlal Results 


Ariel VI 

Balloon 

0.66 ± 0.22 

1.21 ± 0.37 

1.00 ± 0.23 

1.00 ± 0.32 

0.34 ± 0.17 

0.09 ± 0.27 

1.10 ± 0.23 

0.07 ± 0.27 

0.29 ±0.17 

0.94 ± 0.27 












structure, particularly the presence c' a pealt at charge 56. The 
difference between our results and the Ariel VI findings are significant at 
the 75% level using; the same minimuin x 2 irethod discussed previously. 
If we perform the same analysis with the Ariel VI abundances as was 
done in Figures 4.7 and 4.8, we find that the x Z curves resemble ours in 
both the FIP and no FIP cases, although the value of x 2 at the minimum 
is larger for the Ariel VI data in both cases. 


4.7. Summary 

In the previous sections we have shown that a pure r-process source 
is not consistent with our data either with or without adjustments for 
first ionization potential effects. The source composition requires a mix- 
ture of both the r- and the s-process. If one assumes that the cosmic 
ray abundances are altered by the same first ionization potential effects 
as seen at lower charges, then we are consistent with a solar system mix- 
ture of r- and s-process material. The best fit is obtained with a mixture 
which has 


(r/ s) CRS 
(r/ s) ss 


1 . 50 ^ 


where CRS designates the cosmic ray source and SS the solar system. 
On the other hand, if no first ionization potential effects are considered, 
we are most consistent with an s-process dominated mixture having 

= 0.20^o.i4 

and permitting an almost pure s-process source. 

We cannot choose between the FIP and no FIP cases on the basis of 
the present data. Additional steps beyond the scope of this initial study 
will address such matters as obtaining an absolute normalization to the 


(r/ s) C rs 
( r/s)ss 


iron region. This will help determine if the data points are consistent 
with a FIP plot of the sort shown in Figure 4.4. Although our results 
show that we can choose a mixture of r and s-proeess material which, 
when altered by a function of first ionization potential, matches the 
data, the normalization to iron is necessary to establish consistency 
with the other charge regions. If such can be shown to exist, the 
hypothesis of a solar system type source altered by ionization potential 
biases becomes an even more attractive possibility because of the wide 
charge range over which evidence of its effects are seen. As mentioned 
earlier, results in the charge 26 to 40 region, also from the HEAO Heavy 
Nuclei Experiment (Binns et al. 1982), show that the abundances are 
consistent with a solar system type composition altered by FIP. 

The other possibility, that of almost pure s-proeess material, cannot 
be ruled out either. Results from passive balloon borne detectors indi- 
cated the presence of a substantial amount of r-process material in the 
cosmic radiation, based on the detection of a significant flux of trans- 
bismuth nuclei (Fowler et al. 1977). This, result was in agreement with 
the intuitive association of both cosmic rays and the r-process with 
supernovae because of the extremes involved for both: cosmic ray ener- 
gies and r-process neutron fluxes. However, recent results from this 
detector (Binns et al. 1982) have placed a more stringent upper limit on 
the actinide flux with the observation of only one possible actinide 
(88^Z^100) for some 100 platinum-lead (74^Z^87) nuclei, a result which 
is inconsistent with the earlier balloon measurements . The result is, 
nevertheless, still consistent with a solar system type mixture and does 
not require pure s-process material. However, being an upper limit, it 
does not rule out the possibility tither. 

Resolution of the r- and s-process composition of the cosmic radia- 
tion must await a comprehensive picture over the entire charge range to 
remove the free parameters remaining in our analysis. Although we have 
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been using combinations of solar system r- and s-proeess material in an 
attempt to reproduce our observations, it is possible that the cosmic ray 
r-proeess is not the same as that set. a : n the solar system. Nor for that 
matter can we be sure that the s-process is the same. The ledge- 
precipice structure seen in Figure 4.2 is a function of the cross sections 
of magic number isotopes along the s-process path. However, the rela- 
tive levels of the plateaus, or ledges, is a function of the total neutron 
exposure, which may be different for the cosmic rays than for the solar 
system. The detection of a significant difference in abundances between 
elements on the different plateaus or ledges could be an indication of an 
s-process differing from that of the solar system. The most comprehen- 
sive picture, of course, would result from a measurement of the ultra- 
heavy isotopic abundances. But given the combination of the low fluxes 
of the ultra-heavies and the accuracy with which measurements must be 
made in order to resolve isotopes, this goal will not be realized for some 
time. 


136 - 


Appendix A 

The folio-wing is a list of the time periods used in the analysis. Gaps 
of less than 0.1 days were ignored. 


Start End 


Year 

Day 

Year 

Day 

A t 

79 

288.0 

79 

268.3 

0.3 

79 

268.5 

79 

317.0 

48.5 

79 

310.1 

79 

340.0 

22.0 

79 

341.0 

79 

347.0 

6.0 

79 

348.1 

79 

366-0 

17.9 

00 

2.0 

80 

7.0 

5.0 

00 

8.1 

80 

86.1 

78.0 

80 

87.0 

80 

92.0 

5.0 

80 

93.1 

80 

93.9 

0.8 

80 

94.0 

80 

•100.5 

6.5 

80 

101.0 

80 

108.0 

7.0 

80 

109.1 

80 

131.1 

22.0 

80 

133.0 

80 

210.0 

77.0 

00 

211.0 

00 

228.1 

17.1 

80 

230.0 

80 

262.8 

32.7 

80 

263.1 

80 

300.0 

36.9 

80 

300.1 

80 

355.3 

55.2 
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Appendix B 


The decomposition presented in detail in Chapter 4 for the even charge 
elements in the Sn-Ce region was also done for all of the elements with charge 
greater than 50 in order to have a consistent set of abundances to propagate 
over the whole charge range. The r- and s-process abundances that result are 
listed here. As before, the solar system abundances are from Anders and 
Ebihara (19B2), and the s-process is that of Kappeler at aL. (1982). The method 
is the same as detailed in Chapter 4 with s-oniy isotope abundances directly 
from K&ppeler et al. and p-proeess isotopes ignored. The abundances marked 
with b were indicated as having s-proeess branches by Kappeler et al.. This 
branching was ignored in the assigments made here since the effect on the Sn- 
Ce region is small. 
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Element 

— 

Solar 

System 

B 

process 

r 

process 

Element 

Solar 

System 

— 

s 

process 

r 

process 

48 Cd 

1.89 

0.859 

0.691 

68 Er 

0.253 

0.0414 

0.2073 

49 In 

0.184 

0.0617 

0.114 

89 Tm 

0.0386 

0.00539 

0.0332 

50 Sn 

3.B2 

2.123 

1.52 

70 Yb 

0.243 

0.0664 

0.1764 

51 Sb 

0.352 

0.0442 

0.308 

71 Lu b 

0.0370 

0.00563 

D.0313 

52 Te 

4.91 

0.928 

3.96 

72 HI b 

0.176 

0.0656 

0.1097 

53 I 

0.90 

0.0437 

0.856 

73 Ta 

0.0228 

0.00687 

0.0157 

54 Xe 

4.35 

0.66Q 

3.676 

74 ¥ 

0.137 

0.0474 

0.0895 

• 

55 Cs 

0.372 

0.0426 

0.329 

75 Re 

0.0533 

0.00346 

0.0498 

58 Ba 

4.36 

4.11 

0.286 

76 Os 

0.714 

0.0506 

0.6606 

57 La 

0.448 

0.298 

0.150 

77 Ir b 

0.660 

0.0662 

0.594 

56 Ce 

1.18 

0.765 

0.390 

78 Pt 

1.37 

0.0581 

1.312 

59 Pr 

0.174 

0.0764 

0.098 

79 An 

0.186 

0.00783 

0.178 

80 Nd 

0.835 

0.369 

0.390 

80 Hg 

0.52 

0.2068 

0.270 

62 Sm b 

0.262 

0.0608 

0.1810 

81 Tl 

0.184 

0.0831 

0.1011 

63 Eu b 

0.0972 

0.00398 

0.0932 

82 Pb 

3.13 

0.819 

2.310 

84 Gd 

0.331 

0.0395 

0.2892 

83 Bi 

0.144 

0.0183 

0.128 

85 Tb 

0.0589 

0.00323 

0.0557 

90 Th 

0.0420 

0.0 

0.0420 

86 Dy 

0.398 

0.0631 

0.3328 

92 U 

0.0238 

0.0 

0.0238 

87 Ho 


0.0875 

0.00493 

0.0826 
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Element 

A 

Solar 

System 

s- 

process 

process 

4B Cd 

106 

0.0199 

0.0 

0.0 


108 

0.0142 

0.0 

0.0 


110 

0.199 

0.194 

0.0 


111 

0.204 

0.07B9 

0.125 


112 

0.383 

0.207 

0.176 


113 

0.194 

0.0842 

0.110 


114 

0.456 

0.295 

0.161 


116 

0,119 

0.0 

0.119 

49 In 

113 

0.0079 

0.0 

0.0 


115 

0.176 

0.0617 

0.114 

50 Sa 

112 

0.0386 

0.0 

0.0 


114 

0.0256 

0.0 

0.0 


115 

0.0145 

0.0 

0.0 


116 

0.565 

0.458 

0.0 


117 

0.296 

0.105 

0.191 


118 

0.929 

0.655 

0.274 


119 

0.329 

0.156 

0.173 


120 

1.24 

0.749 

0.49 


122 

0.174 

0.0 

0.174 


124 

0.21.5 

0.0 

0.215 

51 Sb 

121 

0.202 

0.0442 

0.158 


123 

0.150 

0.0 

0.150 

52 Te 

120 

0.0045 

0.0 

0.0 


122 

0.123 

0.121 

0.0 


123 

0044 

0.0401 

0.0 


c.\ -JBf , ■ > . . . 

OF PWR vJL'aLi 



134 

0.226 

0.211 

G.O 


125 

0.344 

0.0823 

0.262 


128 

0.918 

C.474 

0.444 


130 

1.56 

0.0 

1.56 


130 

1.69 

0.0 

1.69 

53 J 

127 

0.90 

0.0437 

0.856 

54 Xe 

124 

0.00496 

0.0 

0.0 


128 

0.00483 

0.0 

0.0 


128 

0.0939 

0.108 

0.0 


129 

1.20 

0.0603 

1.14 


13Q 

0.169 

0.176 

n n 

K/.US 


131 

0.941 

0.0619 

0.879 


132 

1.15 

0.254 

0.896 


134 

0.421 

0.0 

0.421 


136 

0.34 

0.0 

0.34 

55 Cs 

133 

0.372 

0.0428 

0.329 

56 @a 

130 

0.00482 

0.0 

0.0 


132 

0.00440 

0.0 

0.0 


134 

0.106 

0.131 

0.0 


135 

0.287 

0.0627 

0.224 


136 

0.342 

0.359 

0.0 


137 

0.488 

0.446 

0.042 


138 

3.13 

3.11 

0.02 

57 La 

138 

0.00040 

0.0 

0.0 


139 

0.448 

0.29B 

0.150 

58 Ce 

136 

0.0022 

0.0 

0.0 


138 

0.0029 

0.0 

0.0 
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140 

1.026 

0.765 

0.261 


142 

0.129 

0.0 

0.129 

59 Pr 

141 

0.174 

0.0764 

0.093 

60 Nd 

142 

0.227 

0.151 

0.0 


143 

0.101 

0.0292 

0.072 


144 

0.199 

0.112 

0.007 


145 

0.0694 

0.0156 

0.0533 


146 

0.144 

0.0607 

0.083 


140 

0.0477 

0.0 

0.0477 


150 

0.0466 

0.0 

0.0466 

flO Cm 

Vlrf L»tu 

1 A A 
i-r-r 

n r»nai 

u ,uuui 

PI p» 

VI<V 

0.0 


147 b 

0.0406 

0.00605 

0.0346 


146 

0.0295 

0.0247 

0.0 


149 

0.0363 

0.00263 

0.0337 


150 

0.0193 

0.0118 

0.0 


152 b 

0.0694 

0.0156 

0.0536 


154 

0.0569 

0.0 

0.0589 

63 Eu 

151 b 

0.0466 

0.00149 

0.0451 


153 b 

0.0506 

0.00249 

0.0481 

64 Gd 

152 

0.00066 

0.0 

0.0 


154 

0.00695 

0.00524 

0.0 


155 

0.0490 

0.00246 

0.0465 


158 

0.0682 

0.0119 

0.0563 


157 

0.0520 

0.00452 

0.0475 


156 

0.0821 

0.0154 

0.0667 


1B0 

0.0722 

0.0 

0.0722 

65 Tb 

159 

0.0589 

0.00323 

0.0557 
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66 Dy 

156 

0.000227 

0.0 

0.0 


156 

0.000398 

0.0 

0.0 


160 

0.00915 

0.00818 

0.0 


161 

0.0756 

0.00232 

0.0733 


162 

0.101 

0.0137 

0.0873 


163 

0.0991 

0.00402 

0.0951 


164 

0.113 

0.0349 

0.0771 

67 Ho 

165 

0.0675 

0.00493 

0.0826 

66 Er 

162 

0.000354 

0.0 

0.0 


154 

0.00395 

n n 

\>.t/ 

G.G 


166 

0.0645 

0.0120 

0.0725 


167 

0.0579 

0.00431 

0.0536 


166 

0.0686 

0.0251 

0.0435 


170 

0.0377 

0.0 

0.0377 

69 Tm 

169 

0.0306 

0.00539 

0 0332 

70 Yb 

186 

0.000320 

0.0 

0.0 


170 

0.00753 

0.00762 

0.0 


171 

0,0350 

0.00415 

0.0309 


172 

0.0532 

0.0145 

0.0387 


173 

0.0394 

0.00686 

0.0325 


174 

0.0768 

0.0331 

0.0437 


176 

0.0306 

0.0 

0.0306 

71 Lu 

175 

0.0359 

0.00457 

0.0313 


176 b 

0.00106 

0.00106 

0.0 

72 Hf 

174 

0.00028 

0.0 

0.0 


176 b 

0.00902 

0.00902 

0.3 
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177 

0.0327 

0 . 0D3B3 

0.0289 


178 

0.0477 

0.0172 

0.0305 


170 

0.0241 

0,00420 

0.0199 


180 

0.0620 

0.0318 

0.0304 

73 Ta 

i8G 

2 . 78e-06 

0.0 

0.0 


181 

0,0226 

0.00687 

0.0157 

74 ¥ 

160 

0.000178 

0.0 

0.0 


182 

0.0360 

0.0169 

0.0191 


163 

0.0196 

0.00980 

0.0098 


184 

0.0421 

0.0207 

0.0214 


106 

0.0392 

0.0 

0.0392 

75 Re 

185 

0.0190 

0.00346 

0.0155 


187 

0.0343 

0.0 

0.0343 

78 Os 

184 

0.000129 

0.0 

0.0 


186 

0.0115 

0.0116 

0.0 


187 

0.0089 

0.00566 

0.0 


188 

0.0954 

0.0128 

0.0826 


189 

0.115 

0.00336 

0.112 


190 

0.189 

0.0172 

0.172 


192 

0.294 

0.0 

0.294 

77 Ir 

191 

0.246 

0.00383 

0.242 


193 b 

0.414 

0.0624 

0.352 

78 Pt 

190 

0.000178 

0.0 

0.0 


192 

0.0107 

0.0106 

0.0 


194 

0.451 

0.0128 

0.438 


195 

0.463 

0,00474 

0.458 


196 

0.347 

0.0300 

0.317 
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190 

0.09B6 

0.0 

0.0966 

79 All 

197 

0.166 

0.00783 

0.176 

60 rig 

196 

0.00076 

0.0 

0.0 


108 

0.052 

0.0103 

0.0 


199 

0.0674 

0.0129 

0.0745 


200 

0.120 

0.0635 

0.056 


201 

0.0686 

0.0329 

0.0357 


202 

0.155 

0.0B72 

0.068 


204 

0.0359 

0.0 

0.0359 

61 T1 

203 

0.0542 

0.0255 

0.0287 


205 

0.130 

0.0578 

0.0724 

62 Pb 

204 

0.0612 

0.0606 

0.0 


206 

0.594 

0.160 

0.414 


207 

0.644 

0.176 

0.466 


200 

1.630 

0.402 

1.426 

63 Bl 

209 

0.144 

0.0163 

0.128 

90 Th 

232 

0.0420 

0.0 

0.0420 

92 U 

235 

0.00573 

0.0 

0.00573 


238 

0.0161 

0.0 

0.0161 
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